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CHAPTER  1 


INTRODUCTION  AMD  LITERATURE  REVIEW 

1.1  Introduction 

Great  strides  have  recently  been  Bade  toward  an  increased  physical 
knowledge  of  turbulence,  especially  in  the  area  of  turbulent  structure.  The 
older,  stochastic  approach  to  the  description  of  turbulent  boundary  layers  is 
being  Modified  as  auch  new  evidence  is  found  to  indicate  the  presence  of  a 
coaplex  deterministic  hierarchy  of  structures.  Soae  basic  structures  have 
been  identified  in  a  variety  of  shear  flows  and  t&ey  appear  to  evolve  through 
tiae  and  space  in  a  quasi- repeat able  manner,  (see,  for  example,  Acarlar  and 
Saith  [1984]).  The  birth  of  the  structures  has  been  linked  to  the  production 
of  turbulent  energy,  their  movement  through  the  boundary  layer  has  been 
associated  with  the  diffusion  of  turbulent  energy,  and  their  breakup  has  been 
linked  with  turbulent  energy  dissipation.  The  importance  of  these  structures 
cannot  be  overemphasized;  to  gain  a  fuller  understanding  of  them,  is  essential 
to  gain  a  fuller  understanding  of  turbulence.  \  v 

Despite  these  observations,  the  actual  form  the  structures  take,  and  how 
they  are  manifested  through  the  boundary  layer,  is  still  a  controversial 
subject.  One  of  the  aims  of  the  present  study  is  to  describe  experimentally 
the  dynamic  behavior  of  these  structures  more  clearly. 

A  further  aim  of  this  study  is  to  exsmine  the  effect  of  compressibility. 
Morkovin  [1962]  hypothesized  that  the  essential  dynamics  of  equilibrium 
compressible  turbulent  boundary  layers  followed  the  incompressible  pattern 
closely,  as  long  as  the  fluctuating  Mach  number  remained  small.  While  this 
hypothesis  has  been  widely  accepted  in  describing  the  t ime-averaged  behavior. 
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the  effect  of  compressibility  on  the  organised  structures  has  not  been  exten¬ 
sively  investigated.  Do  the  structures  appear  the  same,  and  behave  in  the 
ease  way  in  compressible  flow  as  they  do  in  incompressible  flow?  Since  nearly 
all  the  research  on  coherent  structures  has  been  performed  in  incompres¬ 
sible  flows,  we  intend  to  extend  our  knowledge  of  organised  structures  to 
supersonic,  compressible  flows. 

The  work  was  performed  in  a  sero  pressure  gradient  boundary  layer.  As  a 
preliminary,  a  complete  set  of  mean  flow  measurements  was  taken  to  characterise 
the  boundary  layer  used  in  the  study.  The  instantaneous  nature  of  the  turbu¬ 
lence  behavior  was  then  investigated,  and  the  primary  method  of  structure 
eduction  used  conditional  sampling  of  instantaneously  measured  flow  parameters. 
Turbulent  signals  ware  recorded  at  aeveral  locations  simultaneously,  providing 
spatial  and  temporal  information  about  the  structures.  Prior  to  the  condi¬ 
tioning  of  the  signals,  a  complete  time  series  analysis  was  performed,  in¬ 
cluding  cross-correlations  and  cross-spectra. 

The  remainder  of  this  chapter  discusses  the  literature  pertinent  to  this 
experimental  program.  The  equipment,  instrumentation,  experimental  techniques 
and  data  reduction  methods  are  described  in  Chapter  2.  An  overview  of  the 
boundary  layer,  including  mean  flow  aeasuruments  and  the  basic  fluctuating 
quantities,  is  given  in  Chapter  3.  The  two-point  correlation  and  spectrum 
results  are  discussed  in  Chapter  4  and  the  conditional  sampling  results  are 
discussed  in  Chapter  5.  Chapter  6  presents  the  final  discussioo  and  conclu¬ 
sions.  A  tabulation  of  the  mean  flow  data  is  given  in  the  Appendix. 
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1.2  Literature  Review 

The  review  articles  by  Willmarth  [1975],  Cantwell  [1981],  and  Wallace 
[1982]  give  an  excellent  account  of  the  recent  work  on  organized  structures  in 
turbulence.  Willmarth' s  review  deals  with  the  staple  case  of  incompressible 
flow  over  a  smooth,  plane  surface  with  zero  pressure  gradient.  He  considers, 
aeong  other  topics,  spece-time  correlation  measurements,  probe  measurements  in 
the  intermittent  region  and  the  viscous  sublayer,  and  statistical  properties 
of  the  turbulence  field.  Cantwell  and  Wallace  both  concentrate  on  the  orga¬ 
nized  motions  and  flow  structure  present  in  turbulent  fields.  Cantwell 
deals  with  boundary  layers  and  free  shear  flows,  while  Wallace  studies  bounded 
flows  (both  fully  developed  and  transitional).  They  both  give  excellent 
summaries  of  the  work  to  date  along  with  their  own  conceptual  models  for 
organized  motions. 

The  "modern"  history  of  organized  structures  in  turbulent  flow  begins 
with  the  bursting  cycle  first  proposed  in  1967  by  Kline  and  his  colleagues  at 
Stanford  (Kline  et  al.  1967).  A  hydrogen-bubble  wire,  placed  at  various 
positions  near  the  wall,  showed  that  fluid  in  the  viscous  sublayer  did  not 
follow  straight  trajectories;  rather,  it  accumulated  into  alternating  high- 
and  low-speed  areas,  called  streaks.  They  observed  that  the  streaks  had  a 
finite  lifetime  during  which  a  cycle  of  events  was  observed,  called  the 
bursting  cycle.  This  process  is  seen  in  Figure  1.  In  more  detail,  it  consists 
of: 


1)  Accumulation:  The  formation  of  low-  and  high-speed  streaks 
with  a  uniform  spanwise  spacing  at  a  y*  of  3-5. 
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2) 

3) 

4) 

5) 

Kline  el  al.  also  found  that  the  pressure  gradient  had  an  effect  upon  this 
cycle:  a  favorable  gradient  reduced  the  rate  of  bursting,  and  an  adverse 
gradient  increased  the  rate  and  intensity  of  the  cycle.  It  was  suggested 
that  the  pressure  gradient  either  assists  or  hinders  the  lift-up  of  the  streak 
by  effecting  the  inflection  in  the  velocity  profile. 

Kline  et  al.  conjectured  that  the  bursting  phenomenon  plays  a  leading  role 
in  the  production  of  turbulent  energy  and  that  it  controls  the  diffusion  of 
this  energy  from  the  inner  layer  to  the  outer  layer.  Later  research  supported 
this- claim,  most  notably  the  work  of  Corino  and  Brodkey  [1969],  Kim  et  al. 
[1971],  and  Willmarth  and  Lu  [1972].  Corino  and  Brodkey  studied  turbulent 
pipe  flow  and  found  a  mechanism  similar  to  bursting  which  they  termed  "ejec¬ 
tions".  They  noticed  that  these  ejections  occurred  at  a  y*  of  approximately 
10,  and  that  a  "sweep"  followed  the  ejection.  The  sweep  cleans  away  the 
fluid  from  the  previous  ejection  (or  burst).  They  also  found  that  the  inten¬ 
sity  and  number  of  ejections  per  burst  increases  with  Reynolds  number  (a  range 


Wall  Migration:  The  gradual  migration  of  the  low-speed  streaks 
further  out  into  the  wall  layer. 

Lift-up:  The  dramatic  lift-up  of  the  low-speed  streaks  until 
they  penetrate  the  inner  portion  of  the  log-law  region,  thereby 
creating  an  unstable  inflectional  velocity  profile. 

Oscillation:  The  violent,  three-dimensional  oscillation  of 
the  lifted  streak  for  several  cycles. 

Breakup:  The  breakup  of  the  streak  into  much  finer  scales, 
with  a  broad  frequency  content,  which  then  becomes  chaotic. 
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from  2,500  to  50,000).  Most  importantly,  they  observed,  as  Kline  had  sug¬ 
gested,  that  the  bursts  account  for  approximately  70 X  of  the  Reynolds  stress 
level  in  the  boundary  layer  near  the  wall.  Further  work  by  Kia  et.  al.  [1971] 
showed  that  virtually  all  of  the  net  production  of  turbulent  energy  in  the 
range  0  <  y*  <  100  occurs  during  bursts.  Wi 1 Inart h  and  Lu  [1972]  found  that 
significant  contributions  to  the  Reynolds  stress  also  occurred  during  the 
sweep  phase. 

These  early  studies  indicated  that  the  bursting  process  was  essentially  a 
wall  phenomenon .  Further  research,  however,  inferred  a  stronger  level  of 
interaction  between  the  outer  layer  and  the  wall  layer  than  a  random  diffusion 
of  energy  from  inner  to  outer  layer.  Rao  et  al.  [1971]  and  Grass  [1971]  showed 
that  the  mean  bursting  period  scales  with  outer- layer  parameters  (U^  and  5), 
while  the  spanwise  spacing  scales  with  inner-layer  parameters  (u*  and  v).  A 
complex  inter-relationship  between  the  inner  and  outer  layers  is  suggested,  not 
confined  to  a  net  flow  of  fluid  during  the  bursts.  Rao  et  al.  suggested  that 
larger  eddies  from  the  outer  layer  move  down  in  the  boundary  layer  and  "scour" 
the  slow-moving  inner  layer,  creating  regions  of  intense  shear.  The  shear 
excites  and  enlarges  local  instabilities,  thus  triggering  the  burst-sweep 
cycle.  The  cycle  thus  becomes  an  overall  flow  phenomenon  and  can  scale  on  both 
inner  and  outer  flow  variables.  Offen  and  Kline  [1974,1975]  also  showed  that 
the  bursting  sequence  is  affected  by  large-scale  motions.  These  motions  appear 
to  stimulate  a  sequence  of  low-speed  streak  lift-up,  oscillation,  and  breakup. 
Thus  from  beginning  to  end  the  bursting  process  involves  the  interaction 
between  inner  and  outer  layers  (see  Figure  2). 
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The  burst-sweep  cycle  wee  Modelled  by  Many  researchers  us ins  hairpin 
vortices.  These  vortices  take  on  Many  forns,  but  are  basically  A -shaped 
loops  foraed  froa  a  vortex  tube.  In  fact,  Kline  derived  his  bursting  aodel 
froa  studies  of  vortex  loop  structures  developing  in  turbulent  boundary 
layers.  In  these  aodels,  the  aechanisa  of  vortex  stretching  produces  the 
locally  intense  shear  layers  which  then  cause  the  oscillation  and  breakup  of 
the  low-speed  streaks.  The  vortex  stretching,  caused  by  the  interaction  of 
larger  eddies  froa  the  outer  layer  with  the  lifted  low-speed  streaks.  Manifests 
itself  as  a  region  of  concentrated  vorticity  just  outside  the  sublayer.  These 
stretched  vortex  eleaents  undergo  a  rapid  breakup  due  to  their  high  insta¬ 
bility.  The  Models  further  view  the  fomation  of  the  low-speed  streaks  as 
being  due  to  the  concentration  of  low-speed  fluid  between  the  legs  of  the 
vortex.  Two  excellent  exaaples  of  hairpin  vortex  Models  for  near-wall  behavior 
(bursting)  are  provided  by  Hinze  (1975]  and  SMith  [1984]. 

Hinze’s  aodel  begins  with  a  horseahoe  vortex  (a  horseshoe  aay  be  consi¬ 
dered  as  a  short  hairpin  vortex  with  thicker  legs)  which  is  stretched  and 
elongated  by  the  wall  velocity  gradient.  In  accordance  with  the  Biot-Savart 
law,  the  head  of  the  vortex  loop  is  propelled  into  higher  speed  regions  because 
of  self-induction.  This  action  stretches  the  vortex  further,  increasing  the 
vorticity  concentration  in  the  hairpin.  Meanwhile,  low-speed  fluid  is  forced 
between  the  counter-rotating  legs  of  the  vortex,  producing  a  concentrated  shear 
layer.  The  fluid  aurrounding  the  tip  of  the  hairpin  is  subjected  to  the  local 
inflectional  instability  caused  by  this  intense  shear  layer.  The  tip  then 
breaks  down  and  produces  a  turbulent  burst.  The  highly  energetic  turbulent 
fluid  which  results  froa  the  burst  is  convected  downs t rose  and  diffused 
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outward,  thereby  increasing  in  scale  and  decreasing  in  frequency  content.  The 
pressure  waves  associated  with  the  burst  propagate  through  the  boundary  layer, 
setting  up  a  transient  local  pressure  gradient  which  aids  the  aoveaent  of 
"fresh"  fluid  towards  the  wall  -  this  is  the  "sweep".  Hinze  suggests  that  the 
horseshoe  vortices  are  generated  during  the  breakdown  of  the  intense  local 
shear  layers  created  at  the  interface  between  the  inner  and  outer  layers. 

Seith’s  eodel  is  the  result  of  extensive  flow  visualization  studies  using 
side-,  top-,  and  end-views  of  hydrogen  bubble  lines.  To  eodel  the  boundary 
layer  structures  on  a  large  scale,  he  created  vortex  loops  by  using  a  hemi¬ 
sphere  protruding  froe  the  wall  of  the  boundary  layer.  The  heaisphere  served 
as  a  trip  for  the  vortex  sheet  to  roll-up  and  fore  a  horseshoe  vortex  loop; 
see  Figure  3.  When  seen  end-on  (looking  upstreaa),  a  "aushrooaing"  effect  was 
discovered,  which  indicates  a  spreading  of  fluid  in  the  spanwise  direction  and 
infers  that  the  loop  is  leaning  forward.  When  viewed  fraa  the  top,  a  double¬ 
loop  structure  was  observed.  These  double- loop  and  aushrooaing  patterns  are 
the  result  of  a  connected  pair  of  counter-rotating,  streanwise- leaning  vortical 
structures,  that  is,  a  horseshoe  vortex.  While  observing  these  structures, 
Saith  noticed  that  a  horseshoe  vortex  was  always  present  when  a  low-speed 
atreak  evolved,  suggesting  that  the  vortices  are  either  necessary  for  a  burst, 
or  the  result  of  one. 

A  s unwary  of  Saith's  synthesized  aodel  starts  by  considering  a  low-speed 
streak  region.  The  streak  grows  until  the  passage  of  a  disturbance  of  suffi¬ 
cient  size  and  strength  (possibly  caused  by  a  large  eddy  from  the  outer  layer) 
iapresses  a  local  adverse  pressure  gradient  upon  a  portion  of  the  streak.  The 
resultant  local  deceleration  creates  a  three-dieensional  inflectional  profile 
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at  the  interface  between  the  low-speed  streak  and  the  faster  moving  fluid  in 
the  wall  region.  Once  this  inflection  develops,  the  streak  is  extremely 
susceptible  to  local  disturbances  and  instabilities  can  grow.  The  instabil¬ 
ities  cause  the  onset  of  three-dimensional  oscillations  in  the  inflection 
region  which  then  propagate  to  the  vorticity  sheet  encompassing  the  low-speed 
streak.  This  leads  to  the  roll-up  of  the  sheet  into  characteristic  horseshoe 
vortices,  as  we  have  seen  in  Figure  3.  This  full  process  is  shown  in  Figure 
4a  and  4b. 

Once  the  formation  of  the  three-diaensional  vortex  loops  has  begun, 
Biot-Sevart  interactions  between  various  parts  of  the  loop  cause  a  self-induced 
movement  of  the  hairpin  away  from  the  wall.  This  results  in  the  stretching  of 
the  vortex  in  the  streaawise  direction,  which  reinforces  hairpin  formation 
farther  away  from  the  head  of  the  vortex  loop.  A  pair  of  counter-rotating 
legs  with  a  st Teamwise  orientation  is  created,  as  seer  in  Figure  4c.  These 
legs  remain  close  to  the  wall,  but  are  stretched  by  the  strong  velocity 
gradient  there.  Two  phenomena  result:  1)  elevated  viscous  dissipation  due  to 
the  strongly  amplified  vorticity  concentrations  caused  by  the  stretching,  and 
2)  the  creation  of  a  lateral  pressure  gradient,  causing  accumulation  of  low 
speed  fluid  between  the  legs,  which  then  acts  to  perpetuate  the  low-speed 
streaks  (serving  to  close  the  cycle) . 

As  the  head  of  the  hairpin  vortex  is  convected  downstream,  it  creates  a 
strong  local  streamwise  pressure  gradient,  which  causes  the  rapid  ejection  of 
low-momentum  fluid  from  the  low-speed  streak,  as  seen  in  the  end  view  of 
Figure  4c  and  4d.  As  the  burst  event  continues,  the  Biot-Ssvart  effects  cause 
further  distortion  of  the  head  of  the  vortex  and  the  hairpin  aligns  itself  in 


two  characteristic  ways.  First,  the  heads  of  the  hairpin  vortices,  from  those 
just  forming  to  those  near  Mdeatb  ”  due  to  viscous  dissipation,  line  up  at  an 
angle  of  about  15-30°  froa  the  wall.  Second,  the  inclination  of  an  individual 
hairpin  vortex  becomes  approximately  45° ,  the  angle  of  the  principal  axis  of 
strain.  These  two  effects  have  been  seen  in  many  experimental  investigations 
and  are  shown  in  Figure  4d.  Taken  to  the  limit,  the  hairpin  vortex  dissipates, 
leaving  only  the  low-speed  fluid  near  the  wall  (once  the  vortex  loop's  legs). 
Perry  and  Chong  [1982]  point  out  that  "eddy  death"  occurs  because  viscosity 
always  "wins  out”  when  the  structures  are  stretched  at  less  than  exponentially 
increasing  rates. 

Snith’s  aodel  is  probably  the  most  complete  description  of  near-wall, 
organized  structure  behavior.  It  is  a  full  cycle  which  is  able  to  account  for 
much  of  the  experimental  evidence.  Contained  in  the  model  are  Kline’s  bursts, 
Corino  and  Brodkey’s  ejections  and  sweeps,  and  the  typical  shape  of  near-wall 
structures.  Furthermore,  the  model  accounts  for  turbulent  energy  production, 
diffusion,  and  dissipation. 

Researchers  from  Kline  to  Smith  have  concentrated  their  efforts  on  the 
near-wall  organized  structure.  Recently,  however,  research  has  indicated  that 
hairpin  vortices  similar  to  those  connected  with  the  bursting  process  may 
extend  throughout  the  major  portion  of  the  boundary  layer.  Head  and  Bandyo- 
padhyay  [1981]  used  flow  visualization  to  study  a  zero  pressure  gradient 
boundary  layer  over  a  wide  range  of  Reynolds  numbers  (500  <  Reg  <  17,500)  and 
found  significant  Reynolds  number  effects  on  the  boundary  layer  structure. 

The  essential  feature  which  they  discovered  is  the  presence  of  many  hairpin 
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vortices  throughout  the  entire  boundary  layer,  a a  if  the  layer  consisted 
exclusively  of  hairpins  "attached**  to  the  wall  region. 

The  effect  of  the  Reynolds  number  upon  the  character  of  the  vortex  loops 
which  Head  and  Bandyopadhyay  observed  is  seen  in  Figure  5.  These  vortex  loops 
have  undergone  varying  amounts  of  stretching  and  elongation,  depending  on  the 
Reynolds  number.  Apart  from  the  shape  of  the  vortex  loops,  they  found  that 
the  character  of  the  entire  boundary  layer  waa  Reynolds  number  dependent.  At 
higher  Reynolds  numbers  (Reg  >  2000),  the  elongated  hairpin  vortices  originate 
in  the  wall  region  and  extend  throughout  the  boundary  layer.  Individually, 
they  are  inclined  to  the  wall  at  a  characteristic  "eddy  angle"  of  40-50° . 

Larger  notions  appear  to  consist  of  arrays  of  these  hairpins,  with  their  tips 
lying  along  a  straight  line  making  a  smaller  angle  with  the  horizontal  (approx¬ 
imately  20°),  as  shown  in  Figure  6.  This  larger-scale  behavior  was  only 
occasionally  observed  by  Head  and  Bandyopadhyay,  but  it  suggests  that,  at 
times,  the  hairpins  arise  from  the  surface  in  a  systematic  fashion.  They 
speculated  that  as  one  hairpin  leaves  the  surface  it  creates  conditions  which 
are  favorable  for  the  creation  of  another.  At  lower  Reynolds  numbers  (Reg  < 
800)  the  vortex  loops  are  much  less  elongated  and  the  larger  motions  consist 
only  of  isolated  vortex  loops  or  of  several  loops  interacting,  not  of  large 
arrays  of  hairpins  as  found  in  higher  Reynolds  number  flows. 

In  summary,  Head  and  Bandyopadhyay  introduced  a  new  type  of  organized 
structure.  Instead  of  merely  controlling  the  production  and  diffusion  of 
turbulent  energy  (as  with  bursting),  their  proposed  large-scale  structure 
appears  to  describe  the  whole  of  the  boundary  layer.  In  many  ways,  however, 


the  behavior  and  characteristica  of  their  atructure  are  similar  to  the  struc¬ 
tures  used  by  Hinze,  Saith,  and  others  to  aodel  the  bursting  phenomenon. 

Perry  and  Chong  [1982]  have  attempted  to  aodel  the  behavior  of  the  full 
range  of  organized  structures  in  the  boundary  layer  by  eaphasizing  the  theory 
of  vortex  dynaaics  and  their  interactions.  They  proposed  several  shapes  and 
sizes  as  well  as  different  systeaatic  alignaents  of  the  vortices,  ("hierar¬ 
chies").  A  aodel  was  considered  successful  when  it  predicted  the  law  of  the 
wall,  accounted  for  birth  and  death  of  the  eddies ,<•  and  described  the  general 
behavior  of  a  turbulent  boundary  layer  (proper  p.d. f. ,  power  spectra,  etc.). 

The  first  aodel  Perry  and  Chong  put  forth  was  a  simple  aodel  using 
straight  rods  of  vorticity  forming  a  A-vortex.  The  flow  was  modelled  by  a 
random  sampling  of  these  eddies  of  different  sizes,  all  leaning  in  the  down¬ 
stream  direction  at  a  constant  angle  <J>,  (see  Figure  7).  The  model  further 
assumed  that  all  of  the  vorticity  in  the  flow  resided  in  the  vortices,  with  the 
surrounding  fluid  completely  irrotational.  The  law  of  the  wall  can  be  recover¬ 
ed  in  this  instance  by  assuming  conservation  of  circulation  in  the  vortex 
rods,  conservation  of  heat  (in  the  case  of  a  heated  wall),  and  a  "Kline 
scaling"  for  the  lateral  spacing  of  the  vortex  loops  (A  =  100v/uT).  Each 
A-vortex  is  stretched  under  the  action  of  its  surrounding  neighbors,  and  the 
height  of  the  eddy  grows  uniforaly  with  tine.  However,  as  this  stretching 
proceeds,  the  aspect  ratio  of  the  A-vortex  increases,  the  legs  of  the  vortex 
approach  each  other,  and  viscous  diffusion  doainates  over  the  stretching.  The 
resulting  vorticity  cancellation  causes  the  eddy  to  die  at  a  fixed  height 
(which  also  scales  on  v/t^).  However,  it  has  always  been  accepted  that  the 
log-law  region  grows  without  liait  as  61^/v  -*•  88 ,  and  therefore  a  more  advanced 


12 


aodel  is  required  to  correctly  predict  the  outer  layer  behavior. 

A  sore  advanced  aodel  assuaes  geometrically  siailar  hierarchies  of 
eddies,  as  seen  in  Figure  8.  Plane  strain  (froa  neighboring  eddies)  occurs 
within  each  hierarchy  and  the  eddies  grow  froa  their  initial  height  to  a 
height  6,  where  6  is  the  height  of  the  largest  contributing  eddy  in  that 
hierarchy.  In  a  given  hierarchy,  £  is  United  to  the  height  at  which  the 
legs  of  the  A-vortex  begin  to  aerge;  thus,  eddies  larger  than  6  in  a  particu¬ 
lar  hierarchy  are  of  no  consequence,  since  they  are  undergoing  vorticity 
cancellation.  It  is  assuaed  that  each  succeeding  hierarchy  has  a  length  scale 
twice  that  of  its  predecessor.  Since  all  hierarchies  originate  froa  the  sane 
sublayer,  the  characteristic  Mjuap  velocity"  (across  the  original  vortex 

sheet)  is  the  sane  and  scales  with  u  .  The  circulation  thus  doubles  froa  one 

T 

hierarchy  to  the  next,  which  is  a  necessary  condition  for  producing  the  correct 
Reynolds  stress  distribution. 

To  aake  the  aodel  work,  Perry  and  Chong  assuae  that  no  hierarchies,  other 
than  the  first,  have  eddy  heights  less  than  6/2,  where  6  is  the  scale  of  the 
hierarchy  under  consideration.  This  iaplies  that  the  shortest  eddy  in  each 
hierarchy  aust  appear  "out  of  nowhere".  Perry  and  Chong  propose  vortex 
pairing  to  account  for  this  anoaaly;  two  eddies  of  the  largest  height  in  one 
hierarchy  pair  to  give  the  shortest  eddy  for  the  next  hierarchy.  One  aust 
further  assuae  that  half  the  eddies  in  a  given  hierarchy  are  paired,  while  the 
other  half  die  froa  vorticity  cancellation. 

This  aodel  can  be  iaproved  by  using  a  distributed  systea  of  hierarchies, 
as  opposed  to  the  discrete  systea  just  considered.  The  discrete  system 
consists  of  eddies  whose  length  scale  goes  in  gcoaetrical  progression  with  a 
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factor  of  two  froa  one  hierarchy  to  the  next.  This  "quantum  jump"  phenomenon 
can  be  aaoothed  out  by  introducing  randomness  into  the  process  through  jitter 
in  several  quantities.  The  p.d. f.  of  the  eddies  then  consists  of  a  continuous 
distribution  of  scales.  Perry  and  Chong  show  that  either  nodel 
fits  the  law  of  the  wall  and  predicts  other  important  boundary  layer  para¬ 
meters. 

Perry  and  Chong  have  further  discovered,  by  use  of  the  continuous  inverse 
power  law  p.d.f. ,  that  no  wetter  what  shape  the  eddies  have,  and  no  natter  how 
they  are  stretched,  they  will  always  lead  to  a  logarithmic  law,  as  long  as  all 
hierarchies  are  geometrically  similar  and  have  the  same  characteristic  velocity 
scale. 

In  s usury ,  Perry  and  Chong  have  validated,  to  a  certain  extent,  the  work 
of  Head  and  Bandyopadhyay,  Kline,  Smith,  and  others  by  theoretically  predicting 
the  boundary  layer  behavior  using  vortical  structures  similar  to  those  observed 
in  flow  visualization  studies. 

From  this  review  of  previous  work,  it  is  clear  that  our  knowledge  of  the 
existence  and  behavior  of  organized  structures  has  greatly  increased  over  the 
last  twenty  years.  Detailed  measurements  and  visualization  of  these  structures 
have  been  performed  in  the  viscous  layer,  the  wall  layer,  the  log-law  region, 
and  throughout  the  entire  boundary  layer.  Complex  models  have  been  derived  to 
explain  the  structures'  characteristics.  The  wall  layer  has  been  well  docu¬ 
mented  in  regards  to  the  burst -sweep  cycle,  and  the  models  discussed  previously 
are  becoming  widely  accepted.  The  aejor  weakness  of  all  models  is  still  the 
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relationship  between  the  outer-layer  and  the  inner-layer  structures,  and 
further  work  is  obviously  required. 

As  sectioned  previously,  practically  all  research  into  organized  struc¬ 
tures  has  been  confined  to  subsonic  flows.  This  is  not  to  say  that  there  is 
no  evidence  for  coherent  structures  in  supersonic  flows.  For  exaaple,  Head  and 
Bandyopadhyay  [1981]  found  evidence  for  their  hairpin  vortex  hypothesis  froa 
the  schlieren  photographs  of  Decker  and  Neekes  [1976]  and  Decker  [1980].  The 
photographs  show  the  floor  of  a  duct  after  thp  pessage  of  a  shock  wave.  For 
two  different  shock  strengths,  45°  striations  aay  be  found  close  to  the  wall, 
aupporting  the  45°  leaning  vortex  theory.  The  present  authors  have  discovered 
similar  structures  in  two  shadowgraphs  published  by  Van  Dyke  [1982].  Plate  1 
shows  the  turbulent  boundary  layer  on  a  cone  cylinder  Moving  through  the  air  at 
a  Mach  nunber  of  1.84  (photograph  by  A.  C.  Charters,  p.  160);  Plate  2  shows  the 
turbulent  boundary  layer  on  a  body  of  revolution  in  free  flight  at  Mach  2.58, 
(photograph  from  the  U.  S.  Aray  Ballistic  Research  Laboratory,  p.  161). 

Evidence  for  the  45°  turbulent  boundary  layer  structures  which  have  been 
discussed  by  nany  researchers  can  be  clearly  seen  in  these  close-ups. 

Despite  this  flow- visualization  evidence,  the  study  of  turbulent  boundary 
layer  structure  in  supersonic  flow  has  been  largely  neglected.  No  quantitative 
aeasureMents  describing  the  instantaneous  behavior  exist,  and  the  effect  of 
large  density  gradients  and  high  Reynolds  nunber  on  boundary  layer  structure 
are  completely  aisaing.  The  present  study  aims  to  be  the  first  step  in  the 
effort  to  fill  this  gap  in  our  knowledge  of  wall  turbulence. 


Chapter  2 


KXPKBSCirTAL  PBOGRAM:  DfSTBQMNTATION  AND  nCBNIQOKS 

2.1  Overview  of  Program 

To  fulfill  the  alas  of  this  study,  aphasia  was  placed  on  the  siaultaneous 
measurement  of  the  t  lee -dependent  flow  behavior  at  a  nuaber  of  points  in  the 
flow. 

Froa  these  measurements  the  usual  one-point  analyses  can  be  done: 
root-aean-square  values,  probability  density  functions,  power  spectra,  and 
autocorrelation  functions  nay  all  be  coaputed.  In  addition,  two-point  aethods 
such  as  cross-correlations,  coherence  functions,  and  phase  angle  functions  can 
be  used  to  find  length  and  tiae  scale  information  throughout  the  boundary 
layer.  Furthermore,  the  siaultaneous  tiae  histories  can  be  condi tionelly 
saapled  and  examined  for  the  presence  of  characteristic  patterns  in  the 
turbulent  flow. 

The  experiaental  program  consisted  of  four  basic  parts: 

1)  Mean-flow  surveys  of  the  flow  field  to  determine 
boundary-layer  characteristics  and  to  verify  the  two- 
dimensionality  of  the  flow.  Surveys  included:  static 
pressure  profiles,  pitot  pressure  profiles,  and  Preston- 
tube  surveys  (to  deduce  skin-friction). 

2)  Siaultaneous,  multiple  wall-pressure  transducer 
measurements  to  determine  spatially  averaged  convection 
velocities  and  typical  pressure  signals  at  the  wall. 
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3)  Simultaneous,  multiple  hot-mire  anemometer  measurements 
to  determine  typical  eddy  shapes,  sisas,  and  angles  of 
inclination  and  to  determine  typical  ness-f low/velocity 
signatures. 

4)  Simultaneous,  combined  wall-preasure  and  mass-flow 
maaauremants  to  determine  the  connection  between  the 
velocity  field  and  the  pressure  at  tbs  wall. 

2.2  Teat  Bavlronmant 
2.2.1  Wind  Tunnel  Facility 

The  experimental  studies  were  performed  in  the  Princeton  University 
8-inch  by  8-inch  high-Beyno Ids -number,  supersonic,  blowdown  wind  tunnel,  k 
detailed  description  of  this  facility  was  given  by  Ves  and  Bogdonoff  [1971). 

Briefly,  pressurised  air  is  supplied  by  four  Worthington  four-stage 
compressors  (from  atmospheric  to  3200  psia)  and  stored  in  four  tanks  with  a 
total  capacity  of  2000  cubic  feet.  A  hydraulically  controlled  valve  releases 
the  pressurised  air  into  the  settling  chamber  at  a  stagnation  preesure  preset 
between  60  and  600  pei.  The  settling  chamber  temperature  decreases  by  several 
degrees  Kelvin  during  each  run  due  to  expansion.  Pram  the  stilling  chamber, 
the  air  ia  expanded  through  a  convergent-divergent  nossle  to  a  nominal  Mach 
number  of  2.9  in  the  working  section.  The  working  section  is  composed  of  three 
interchangeable  35.8  inch  long  constant  area  test  sections  with  a  cross-section 
of  8  inch  x  8  inch.  The  air  is  than  exhausted  to  atmoephere  through  a  diffuser 


section. 


2.2.2  That  Section 


Tost#  were  conducted  in  the  floor  boundary  layer  in  the  second  section 
dowostreen  tram  the  nossle.  For  convenience,  a  coordinate  systen  wee  refer¬ 
enced  to  the  start  of  this  section  with  X  the  streaawise  coordinate,  Z  the 
•penwise  coordinate  (with  Z  *  0  on  the  centerline)  and  T  the  coordinate  nornal 
to  the  well.  The  origin  of  this  coordinate  systen  is  68.5  inches  dowostreen 
froa  the  nossle  throat  (see  Figure  9),  neaeured  along  the  tunnel  wall.  To 
■ensure  the  static  pressure  with  a  static  pressure  tap  or  niniature  wall- 
pressure  transducer,  a  two  inch  lastrvwntstion  plug  was  inserted  on  the 
centerline  of  the  tunnel  at  X  *  18.5  inches.  To  survey  the  floor  boundary 
layer,  probes  were  inserted  froa  one  of  two  different  ceilings.  For  asasure- 
aents  taken  along  the  tunnel  centerline,  a  ceiling  with  a  sirs— iilse  slot  was 
used.  For  aeesuronsnts  taken  off  the  centerline,  a  ceiling  was  used  which 
could  access  any  atreanwiae  location,  yet  also  gave  lateral  nobility. 

2.2.3  Teat  Conditions 

For  all  teats  the  stagnation  pressure  was  100  pels  ±  0.5k  and  the  stag¬ 
nation  teapereture  was  noninally  270  Kelvin.  For  hot-wire  aeaaurenents,  the 
running  tiae  was  approximately  two  ainutes,  with  a  stagnation  tenperature  drop 
of  3-4%;  other  runs  were  considerably  shorter  and  the  taaperature  drop  was  not 
“  greet.  The  flow  had  a  typical  freestrean  Mach  maber  of  2.87  (+  0.01)  with 
*  unit  Reynolds  maber  of  6.5l+'.7/a  (+  4%).  The  walls  were  approx iaately 
adiabatic  aid  the  frees treaa  turbulence  level  was  1-1.5%.  All  of  the  instan¬ 
taneous  flow  aeaauraaanta  were  centered  about  X  *  18.5  inches,  Z  *  0  inches. 
Conditions  at  that  point  ore  given  in  Table  1. 


18 


2.3  Data  Acquisition 

2.3.1  Data  Acquisition  Hardware 

Separate  data  acquisition  systems  were  used  for  aean  flow  surveys  and 
for  instantaneous  measurements . 

2.3. 1.1  Preliminary  Naan  Flow  Surveys 

For  the  aean-flow  surveys,  the  data  acquisition  system  consisted 
of  a  Hewlett-Packard  1000  minicomputer,  s  Preston  Scientific  GMAD-4  analog-to- 
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digital  converter,  and  a  Hewlett-Packard  2240  Measurement  and  Control  Pro¬ 
cessor.  The  GMAD-4  is  a  sequential,  multi-channel,  14  bit  (plus  sign)  analog- 
to-digital  converter  with  an  input  range  of  +  10  volts  and  a  maximua  sampling 
frequency  of  50  kHz.  The  Measurement  and  Control  Processor  was  interfaced  with 
the  minicomputer  and  controlled  probe  movement  during  flow  surveys. 

Very  briefly,  the  output  voltage  from  each  transducer  (stagnation  pressure 
and  temperature,  position,  and  quantity  of  interest)  was  digitised  by  the 
GMAD-4,  converted  from  counts  to  real  veluee  by  the  appropriate  calibration 
curve,  and  then  stored  on  the  coaputer  bard  disc. 

2.3. 1.2  Fluctuating  Data  Surveys 

For  the  fluctuating  data  collection,  a  much  faster  system  was  used,  with 
sampling  rates  of  up  to  one  million  hertz  on  each  of  four  simultaneous  chan¬ 
nels.  The  system  used  consisted  of  a  VAX  11/750  minicomputer  and  CAMAC 
(Computer  Automated  Measurement  and  Control)  data  acquisition  system.  The 
CAMAC  system,  obtained  from  LeCroy  Inc.,  consisted  of  a  programmable  amplifier, 
an  analog-to-digital  converter,  and  memory.  All  units  were  designed  to  handle 
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froa  one  to  four  channels  of  data  simultaneously,  and  they  were  installed  in  a 
Kinetic  Systems  Crate.  Interfacing  with  the  VAX  was  accomplished  through  a 
crate  controller. 

The  analog  data  signal  entered  the  CAMAC  system  through  the  8100  Dual 
Programmable  Amplifier,  equipped  with  adjustable  gains  from  0.2  to  100.  The 
D.C.  bias  was  adjusted  to  zero  at  the  specified  gain,  and  the  amplifier 
roll-off  frequency  was  greater  than  1  MHz  on  each  channel.  The  amplified 
signals  were  then  sampled  simultaneously  (less  than  5  nanosecond  uncertainty) 
and  digitized  at  a  rate  of  1  MHz  by  the  TD8210  Waveform  Analyzer,  with  a 
resolution  of  9.8  millivolts/count  (10  bits)  over  an  input  range  of  +  5  volts. 
All  four  of  the  TD8210  channels  had  slightly  different  gains,  which  were 
accounted  for  in  the  data  reduction  programs.  These  digitized  signals  were 
then  temporarily  stored  in  a  96000  word  (98304  data  points)  memory  formed  by 
cascading  three  8800A  Memory  Modules.  When  this  memory  became  full,  the  data 
mas  unloaded  (stripping  one  channel  at  a  time)  to  disc  storage  on  the  VAX  (as 
binary  representations  of  A/D  counts).  The  nearory  modules  allowed  a  maximum, 
continuous  record  length  of  24576  data  points  per  channel  for  each  of  four 
channels. 

Another  integral  part  of  the  high-speed  data  acquisition  system  was  a 
number  of  Ithaco  Model  4213  analog  filters  with  a  roll-off  of  -80  db/decade 
(four  pole  Buttemorth  filters). 

Since  much  of  the  experimental  analysis  required  comparisons  between 
channels,  no  phase  shift  among  the  channels  could  be  tolerated.  Therefore, 
smch  CAMAC  channel  and,  more  importantly,  each  Ithaco  filter  was  checked  to 
ensure  a  zero  relative  phase  shift  in  the  frequency  range  of  operation.  This 
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wu  accomplished  by  inputting  a  sine  wave  to  two  different  channels,  and 
monitoring  the  output  for  any  relative  shift  of  signals.  Bach  of  the  channels 
and  filters  which  were  used  was  found  to  have  virtually  no  relative  phase 
shift. 

2.3.2  Data  Acquisition  Procedures 

2.3.2. 1  Preliminary  Measurements 

Each  nean-flow  data  point  consisted  of  four  measurements :  stagnation 
pressure  and  stagnation  temperature  in  the  settling  chamber,  a  linear  distance 
(fron  the  floor,  for  instance)  and  the  quantity  of  interest  (pitot,  static,  or 
Preston-tube  pressure).  At  the  start  of  the  day,  each  transducer-A/D  converter 
pair  ms  calibrated  by  inputting  a  known  signal  (froai  a  laboratory  standard) 
for  at  least  five  values  in  the  range  of  interest.  Calibration  coefficients 
were  deduced  from  a  least-squares  linear  fit  to  the  data  pairs  of  input  signal 
and  corresponding  sampled  value.  Only  calibrations  with  a  standard  deviation 
of  less  than  1*  were  accepted. 

The  stagnation  pressure  was  measured  with  a  Pace  500  psi  transducer 
referenced  to  atmosphere.  It  was  calibrated  against  a  Heise  0-500  psi  pressure 
gauge,  accurate  to  0.2*  of  its  full  scale  deflection,  in  10  psi  increments 
fron  80  to  120  psi.  Pitot,  static,  and  Preston-tube  pressures  were  measured 
with  10,  50,  or  100  psi  Druck  transducers  referenced  to  vacuum  and  calibrated 
against  a  Wallace  and  Tieraan  standard  gauge.  The  stagnation  temperature  was 
measured  by  a  chroael-alumel  thermocouple  with  an  ice  bath  reference.  The 
calibration  was  performed  against  a  Tine  Electronic  millivolt  source  accurate 
to  0.05*  of  full  scale  (9.99  millivolts).  Linear  distances  were  measured 
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using  potent toasters  calibrated  against  dial  indicator  gauges,  narked  in  .001 
inch  increments. 

All  pressure  probes  were  designed  to  be  long  and  slender  in  order  to 
aininise  local  interference.  The  static  pressure  probes  were  snail,  cone- 
tipped  cylindrical  probes  nade  frost  0.85  —  dianeter  hypodenic  tubing. 

IWo  0.25  ■  orifices  were  located  ten  dia— ters  behind  the  tip,  circunferen- 
tially  opposite  each  other  in  the  horizontal  plane.  The  pitot  pressure  probe 
had  a  0.18  an  flattened  tip  with  a  0.08  ■  opening.  The  Preston  tube  had  a 
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circular  cross-section  with  an  inner-to-outer  dianeter  ratio  of  0.6. 

As  each  probe  was  traversed,  a  single  pressure  neasurenent  was  taken  at 
each  location.  No  averaging  was  done,  but  a  delay  of  300  nillieeconds  was 
used  to  allow  the  pressure  transducer  to  equilibrate.  This  delay  was  found  to 
give  satisfactory  results  under  all  operating  conditions  encountered  in 
this  experiaent. 

2. 3. 2. 2  Fluctuating  Nail-Pressure  Measure— nts 

Measurements  of  the  wall-pressure  fluctuations  were  nade  using  four 
indent! cal  miniature  differential  pressure  transducers  nanufactured  by  Kulite 
Seal conductor  Products  Inc.,  Model  X<XH)62-25-D.  Bach  transducer  had  a  0.71 
—  dia— ter  silicon  sensing  eleaent  on  which  a  fully  active  Wheatstone  bridge 
was  bonded  atonically.  The  transducer  emitted  a  voltage  proportional  to  the 
difference  between  the  pressure  applied  to  the  silicon  diaphre *n  and  a  refer¬ 
ence  pressure.  For  this  experiaent,  the  reference  pressure  was  set  to  the 
local  —an  pressure  (3.34  psi).  Therefore,  the  seen  pressure  was  eliminated 
from  the  signal  and  the  resolution  of  the  fluctuating  coaponent  could  be 
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increased  dramatically.  The  natural  frequency,  aa  quoted  by  the  manufacturer , 
was  500  KHz.  The  best  estimate  for  the  usable  frequency  range  was  approx¬ 
imately  0-80  KHz. 

The  transducers  were  calibrated  statically  every  three  or  four  runs  at 
the  operating  temperature,  by  applying  a  known  pressure  to  the  transducer 
diaphragm  and  monitoring  the  output.  The  calibrations  were  carried  out  at  the 
operating  temperature  because  it  was  found  that  the  resulting  calibration 
coefficients  had  a  strong  dependence  on  the  stagnation  temperature.  Once 
again,  only  calibrations  with  a  standard  deviation  of  leas  than  one  percent 
were  accepted.  It  should  be  noted  that  shock  tube  tests  by  Raman  [1974]  have 
shown  that  transducers  of  this  type  have  dynamic  calibrations  only  a  few 
percent  lower  than  those  obtained  statically. 

Four  pressure  transducers  were  mounted  in-line  in  a  cylindrical  plug  (see 
Fig.  10a),  which  was  then  fitted  in  the  test  section  floor  with  its  center  at 
X  =  18.5  inches,  Z  =  0  inches.  The  separation  distance  between  each  transducer 
was  0.2  inches,  that  is,  transducers  were  located  at  X  =  18.2,  18.4,  18.6,  18.8 
inches.  The  plug  could  be  rotated  through  360°,  yawing  the  transducers 
relative  to  the  flow. 

Hanly  [1975]  has  shown  that  the  flushness  of  the  transducers  is  an 
important  parameter  in  measuring  fluctuating  surface  pressure  accurately. 
Accordingly,  the  transducers  were  adjusted  to  less  than  0.0026  under  the 
floor  surface  using  a  microscope. 

Due  to  the  low  output  voltage  level  of  the  transducers,  each  signal  was 
amplified  by  a  two-stage  amplifier  circuit  for  a  total  gain  of  1000.  These 
amplifiers  were  found  to  have  roll-off  frequencies  of  at  least  400  KHz. 
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The  signals  were  then  band-pass  filtered  by  the  Ithaco  analog  filters.  The 
high-pass  was  set  to  250  Hz  to  reduce  electronic  noise  and  low  frequency 
tunnel  noise.  The  r.a.s.  noise  in  the  systea  (without  the  tunnel  running)  was 
found  to  be  .004  psi,  giving  a  typical  signal-to-noise  ratio  of  12.  The 
low-pass  was  first  set  to  500  KHz,  which  is  one-half  the  sampling  rate,  to 
reduce  aliasing  of  the  signal.  At  such  a  high  cut-off  frequency,  however,  a 
resonance  of  the  Kulite  transducers  at  480  KHz  was  producing  a  very  large  spike 
in  the  energy  distribution  (as  found  froa  the  power  spectra).  Since  the  gain 
of  the  transducers  is  constant  to  only  80  KHz,  it  was  decided  to  low-pass  at 
125  KHz  to  reduce  the  effect  of  the  reaonance.  Even  so,  the  spectra  of  the 
transducers  was  unreliable  down  to  40  KHz.  The  effect  of  this  will  be  seen  in 
Chapter  3.  After  the  signals  were  filtered,  they  were  fed  into  the  CAMAC 
systea  where  they  were  anplified  by  another  factor  of  5,  digitized  at  a  rate  of 

1  MHz,  and  written  to  disc  on  the  VAX.  See  Figure  11  for  a  flowchart  of  the 
pressure  data  acquisition  systea.  Given  a  typical  pressure  fluctuation  of  0.05 
psi,  the  combined  gain  of  5000  results  in  an  A/D  resolution  of  60  counts  for 
that  pressure  peak. 

Data  was  obtained  in  files  of  three  records,  each  record  containing 
24,576  points  per  channel.  Convergence  of  Bean  pressure,  r.a.s.  pressure, 
probability  density  distribution,  and  power  spectrua  was  easily  achieved  under 
such  sampling  procedures. 

Phase  shifts  between  the  transducers  were  determined  indirectly  by  taking 
one  run  with  the  transducers  aligned  with  the  flow,  another  after  rotating  the 
plug  180°,  and  then  comparing  the  cross-correlations.  A  aaxiaua  shift  of 

2  aicro-seconds  was  observed  between  channels;  i.e.,  the  location  of  the 
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maxima  of  the  cross-correlation  differed  by  a  maximum  of  2  micro-seconds  when 
the  plug  was  rotated.  The  implications  of  this  result  are  discussed  further 
in  Chapter  3. 

2. 3. 2. 3  Fluctuating  Mass-Flow  Measurements 

A  DISA  55M10  Constant  Teaperature  Hot-Wire  Anemometer  was  used  to  measure 
the  instantaneous  mass  flux.  The  hot-wire  probes  were  constructed  by  electro¬ 
plating  5  micron  diameter  tungsten  wire  with  copper  and  soft-soldering  the 
wire  onto  the  probe  prongs  (2.5  -  3  me  separation).  The  wires  were  then 
etched  using  a  dilute  sulphuric  acid  solution  to  expose  an  active  portion  of 
tungsten  wire  0.8  -  1.0  am  in  length.  A  small  amount  of  slack  was  introduced 
into  the  active  length  of  the  wire  to  avoid  strain-gauging  at  high  frequencies. 
This  type  of  probe,  first  recoMended  by  Kovasznay  [1950]  and  then  critically 
examined  by  Saits  et  al.  [1983],  gives  reasonable  wire  lifetimes,  minimizes 
interference  from  bow  shocks  emanating  from  the  prongs,  and  yields  a  good 
frequency  response. 

For  multi-wire  runs,  a  special  hot-wire  support  was  designed  to  hold 
four  normal  wires  in  pairs  of  two,  one  above  the  other,  and  the  two  pairs  of 
wires  could  be  moved  relative  to  each  other  vertically  (see  Fig.  10b).  No 
runs  were  actually  made  using  all  four  wires.  However,  two  runs  were  made  with 
a  double  normal-wire  probe:  one  with  a  y-separation  of  2.38  am,  the  other  with 
2.58  wm.  Two  runs  were  also  made  with  a  triple  normal-wire  probe:  the  sepa¬ 
ration  distances  were  2.58  mm,  3.17  am,  and  5.75  mi. 

In  a  constant  temperature  anemometer  system,  the  output  signal  contains 
contributions  from  mass-flow  fluctuations  and  stagnation  temperature  flue- 
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tuations.  The  hot-wire  response  to  both  contributions  varies  with  the  overheat 
ratio,  t,  defined  as: 


T  -  (T„-T.)/T0 

where  T„  is  the  wire  temperature,  To  is  the  flow  stagnation  temperature,  and 
T«  is  the  wire  recovery  teaperature.  The  contribution  due  to  aass-flow 
fluctuations  and  temperature  fluctuations  could  be  separated  by  operating  the 
wire  over  a  range  of  overheat  ratios.  However,  the  constant  temperature  system 
is  inherently  unsuitable  for  operation  at  low  overheat  ratios,  and  it  is 
therefore  unsuitable  for  measuring  stagnation  temperature  fluctuations.  Smits 
et  al.  [1983]  have  shown  that  at  higher  overheat  ratios  (between  1.0  and  1.3), 
the  contribution  of  the  temperature  fluctuations  is  small  enough  so  that  it  can 
be  neglected  and  the  hot-wire  anemometer  is  then  sensitive  only  to  mass-flow 
fluctuations.  In  addition  to  diminishing  the  effect  of  the  stagnation  temper¬ 
ature  fluctuations  upon  the  output  signal,  increasing  the  overheat  ratio 
improves  the  frequency  response  of  the  wire.  The  upper  limit  of  the  overheat 
ratio  is  determined  by  the  oxidation  of  the  tungsten  wire.  With  a  stagnation 
temperature  of  270K,  the  maximum  overheat  ratio  is  approximately  1.0,  and  this 
value  was  used  throughout  the  present  series  of  tests. 

A  small  Mach  3  pilot  tunnel  war  used  to  calibrate  the  wires  and  adjust 
their  frequency  response.  The  frequency  response  was  determined  by  operating 
the  wire  in  the  freestream  of  the  pilot  tunnel  (at  the  prescribed  stagnation 
pressure  for  all  tests,  100  psi)  while  superimposing  a  square-wave  of  small 
amplitude  on  the  wire  voltage.  Perry  and  Morrison  [1971]  have  demonstrated 
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the  usefulness  of  this  aethod  for  subsonic  flow,  while  Bonnet  [1982]  hss 
extended  the  aethod  to  coapressible  flows.  The  frequency  response  is  optimized 
by  adjusting  the  fain  of  the  anenoaater,  the  level  of  the  high  frequency 
filter,  and  the  inductance  of  the  circuit  until  the  transient  response  of  the 
systea  is  ainisized.  An  acceptance  criteria  for  each  wire  was  that  the  upper 
roll-off  frequency  be  at  least  120  KHz. 

The  spatial  resolution  of  the  hot-wire  probe,  as  determined  by  the  active 
length  of  the  wire,  restricts  the  frequencies  which  can  be  discerned.  Since 
the  active  portion  is  larger  than  the  smallest  turbulence  scales,  the  hot  wire 
distorts  the  contribution  of  the  high  frequency  end  of  the  spectrum,  Wyngaard 
[1968]  suggested  that  the  aeasured  value  of  the  one-diaensional  spectra  falls 
to  one-half  of  its  value  at  a  wave  nuaber  of  2.1/t  (where  l  is  the  active 
length  of  the  wire).  Under  test  conditions,  for  s  length  of  0.8  an,  this 
frequency  was  250  KHz;  since  this  is  greater  than  the  frequency  limitation  of 
the  entire  systea,  it  is  not  the  aost  restrictive  factor. 

To  calibrate  the  hot  wires  for  aass-flow  sensitivity,  the  wire  temperature 
(and  therefore  the  overheat  ratio)  was  set  by  selecting  the  wire  operating 
resistance,  and  the  aass-flow  rate  was  varied  by  changing  the  stagnation 
pressure  over  the  expected  operating  range  (in  practice,  fraa  80  to  170  psi). 
The  wall  static  pressure,  stagnation  temperature,  and  anemometer  output  voltage 
were  all  Monitored  along  with  the  stagnation  pressure.  The  aass-flow  rate  was 
calculated  and  the  calibration  data  was  fitted  with  the  following  fora  of 
King’ a  Law: 


B2  =  L  +  M(pu)n 


Fro*  this  relation,  the  calibration  coefficients  L  and  M  could  be  calculated 


along  with  the  aass-flow  exponent,  n,  simply  by  fitting  a  linear  curve  to  E2 
vs.  (pu)n.  The  stagnation  temperature  uaually  varied  by  several  degrees 
during  a  calibration  run,  and  the  coefficients  L  and  M  were  corrected  in  an 
iterative  aanner  (as  described  by  Saits  et  al.  (1983}).  The  calibration 
tesperature  to  which  all  data  points  were  referenced  was  270X.  A  value  for  n 
of  0.55  proved  to  give  excellent  results  for  all  calibration  data.  For  a 
typical  calibration  curve  see  Figure  12. 

Boundary- layer  traverses  were  aade  at  aany  streaawiae  and  spanwise 
stations.  These  traverses  were  aade  using  single  normal  wires  or  multiple 
normal  wires  arranged  in  an  array.  The  voltage  output  froa  each  anesnaeter 
was  separated  into  aean  and  fluctuating  components,  by  aeans  of  the  Ithaco 
analog  filters,  for  increaaed  resolution  of  the  fluctuating  signal.  The  aean 
component  was  obtained  by  low-pass  filtering  at  10  Hz,  digitizing  the  signal 
at  50  KHz  with  the  GMAD/4  analog-to-digital  converter,  and  then  writing  the 
data  to  disc  on  the  VAX.  The  stagnation  teaperature,  stagnation  pressure,  and 
a  linear  measurement  of  the  hot-wire  probe  froa  the  floor  were  similarly 
monitored  by  the  CMAD/4. 

The  fluctuating  coaponent  of  each  aneaoaeter  signal  was  band-pass  filtered 
from  10  Hz  to  500  KHz,  amplified  by  a  factor  of  5  and  sampled  by  the  CAMAC 
system  at  1  MHz  per  channel.  The  low  frequency  limit,  10  Hz,  was  chosen  to 
match  the  mean  signal  filter,  and  the  high  frequency  limit,  500  KHz,  was  chosen 
as  one-half  the  sampling  rate  to  avoid  aliasing  errors  of  the  spectrum.  See 
Figure  13  for  a  flow  chart  of  the  hot-wire  data  acquisition  system.  Data  from 
each  point  in  the  traverse  was  written  to  a  separate  file  containing  one  record 
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of  24,576  points  por  channel.  This  record  length  gave  adequate  convergence  of 
the  r.a.a.  valueo  (within  lk  at  y/6  »  0.5),  the  probability  density  function, 
and  the  power  spectnai.  The  total  instantaneous  values  of  the  nans  flow  were 
found  by  adding  the  seen  voltage  to  the  fluctuating  voltage  and  inverting  the 
calibration  curve  directly,  thereby  avoiding  the  uee  of  the  linearised  sensi¬ 
tivity  coefficients.  Before  inverting  the  calibration  curve,  however,  both  the 
total  voltage  and  the  tine- averaged  naan  voltage  had  to  be  corrected  for 
changes  in  the  stagnation  t separatum  between  the  actual  run  and  the  cali¬ 
bration  according  to 
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where  Tc  is  the  reference  calibration  stagnation  tanperature  (2701)  and  To  is 
the  stagnation  tanperature  during  the  run.  The  instantaneous  values  of  the 
seas- flow  fluctuations  warn  than  obtained  by  subtracting  the  tine-averaged 
aaas  flow  froa  the  total  instantaneous  boss  flow. 


2. 3. 2. 4  Fluctuating  Mass  Flow/Kal  1-Praeeura  Noasurananta 
By  caabining  the  use  of  wall-pressure  transducers  and  hot  wires,  simul¬ 
taneous  ■aaaursnants  of  the  instantaneous  wall  pressure  and  instantaneous  nass 
flow  were  aade.  The  hot  wires  warn  placed  at  different  points  in  the  flow 
(relative  to  the  wall -pres sum  transducers)  to  obtain  s  wide  spatial  resolution 
of  the  flow  field. 

Several  different  configurations  of  pressure  transducers  and  hot  wires 
warn  used: 
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1)  Three  Kulites  aligned  with  the  flaw  (separation  distances  of  0.2  and 
0.4  inches)  end  single  Dorsal  hot-wire  surveys  at  four  upstreaa  stations 
(ranging  to  -1.0  inch  from  the  farthest  upstreaa  Kulite)  and  tan  downs t re sa 
stations  (ranging  to  2.0  inches  downs tress  of  the  farthest  downs t rasa  Kulite). 

2)  Kulitas  aligned  as  above,  but  hot-wire  traverses  at  five  spanwise 
stations  (ranging  to  0.7S  inches  laterally  away  fraa  the  line  of  pressure 
transducers) . 

3)  Two  Kulites  aligned  with  the  flow  and  double  no  real  wire  traverses  at 
several  downs tresa  stations. 

Once  again,  the  aatter  of  relative  phase  shift  between  the  signals  is  an 
iaportant  consideration.  The  best  way  to  determine  the  phase  shift  would 
involve  having  a  step  input  aeasured  simultaneously  by  both  a  hot  wire  and 
pressure  transducer.  By  performing  a  Fourier  transform  on  the  output  signals 
to  determine  frequency  content,  and  examining  this  spectra,  the  limit  of  the 
response  can  be  found.  Alternatively,  by  examining  the  frequency  limitations 
of  both  systems,  this  limit  can  be  found  approximately.  As  we  saw  previously, 
the  frequency  response  of  the  hot-wire  s ye tea  was  always  flat  to  at  least 
120  KHz,  whereas  the  Kulite' s  maximum  frequency  response  was  limited  to 
approximately  90  KHz.  However,  the  spectra  of  the  Kulites  above  40  KHz  was 
found  to  be  unreliable  (this  will  be  seen  in  Chapter  3).  Hence,  the  upper 
limit  on  the  coabined  measurements  is  approximately  40  KHz. 


2.4 
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2.4.1  Tim  Sorias  Analysis 

Once  tho  aeon  and  r.a.s.  were  deterained  for  tha  aass-flow  and  the  wall 
praaaura,  tins  aariao  analysis  was  usad  to  furthar  axaain*  tba  data.  This 
analysis  includad  tha  probability  dansity  function,  tha  powar  spactrun,  and 
for  aulti-channal  analysis  tha  cross -corra 1st ion,  tha  coapactrun,  tha  cobar- 
onca,  and  tho  phasa  angle. 

Before  any  data  analysis  was  parfornad,  sins  waves  and  square  waves  ware 
used  to  det arsine  that  all  routines  wara  operating  correctly. 

Tba  reduced  fluctuating  data  was  analysed  as  1024  point  enseables  in  the 
tins  aariao  program.  Thera  were  several  rasa  one  for  choosing  this  record 
length,  including  the  mount  of  confutation  tine  required  to  calculate  Fourier 
transform  on  tha  anaanbla.  Another  reason  cam  froa  considerations  involving 
the  cross-correlation  function.  Tba  cross-correlation  progrm  was  tested  (on 
all  three  types  of  data:  pres sure/ pressure ,  mss  flow/aass  flow,  and  pressure/ 
aass  flow)  using  record  lengths  of  256,  512,  1024,  2048,  4096,  8192,  and 
24,576  points.  The  aaxiaua  value  of  the  cross-correlation  always  occurred  at 
the  sane  value  of  the  tine  delay  (regardless  of  the  ana cable  length)  yet  the 
peak  value  varied  by  ♦  15k.  Tba  trend  was  for  lower  peak  values  at  amll 
record  lengths,  a  aaxiaua  at  1024  points,  followed  by  lower  peak  values  again 
at  larger  ana sable  lengths.  Ibis  is  probably  caused  by  "phase  jitter",  which 
affectively  smart  the  peak  value  due  to  the  large  nuaber  of  realisations  in 
longer  unstable  lengths.  Since  the  aaxiaua  value  of  the  peak  occurred  at  a 
record  length  of  1024  points,  this  ainiaised  the  phase  jitter  yet  was  long 
enough  to  achieve  a  snoot h,  well-averaged  peak  value.  Also  considered  was  the 
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fact  that  the  langth  of  tha  enaeeble  determines  tha  lowest  fraquancy  for  the 
tiaa  aariaa  analysis.  With  a  record  langth  of  1024  points,  this  lower  Unit  is 
approximately  1  KHz.  For  tha  present  boundary  layer,  vary  United  energy  is 
contained  below  1  KHz,  and  therefore  a  record  length  of  1024  points  was 
considered  to  be  acceptable. 

In  each  tine  series  analysis,  except  for  the  probability  density  function, 
fast  Fourier  transform  were  used  to  iaprove  the  speed  of  the  conputation.  To 
coapute  the  Fourier  coefficients  a  progran  called.  FOURT  (written  by  Noman 
Brenner,  Charles  Rader,  and  Ralph  Alter,  MIT  Lincoln  Laboratory,  August  1967), 
was  used,  which  enploys  the  Cooley- Tukey  fast  Fourier  tr ana fora  net hod. 

2.4. 1.1  Probability  Density  function  (pdf) 

The  probability  density  function  provides  infomation  on  the  anplitude  of 
the  fluctuating  data;  it  describes  the  probability  that  the  data  will  lie 
within  a  defined  range  of  values  at  each  instant  in  tine.  We  define  the  pdf 
as 

p(x)  -  lin 
Ajr*o 

If  the  data  is  noraally  distributed,  the  pdf  will  have  a  Gaussian  distri¬ 
bution.  Hie  Mali  pressure  followed  this  closely,  as  shown  in  Figure  14. 

Effects  such  as  interaittency  produce  a  skewed  pdf,  as  is  evident  near  the 
edge  of  the  boundary  layer  with  the  aass-flow  data  (see  Figure  15). 

In  detemining  the  pdf,  100  cells  (or  classes)  of  data  were  chosen  over 
SlH  standard  deviations,  giving  a  ceil  width  of  0.06o.  The  probability 
^nsity  function  was  then  estiaated  by 


ProbCx<x(t)<x+Ax, 

Ax 
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P<«>  •  O  . 

where  N  is  tbs  total  mafcor  of  data  pointa,  W  ia  tbs  call  width,  and  Mb  is  the 
nusber  of  data  points  which  fall  within  tba  range  X  *  M/2. 

2.4. 1.2  Fewer  Spectral  Density  Function 

Tbs  power  spectral  density  function  gives  infornation  as  to  tbs  frequency 
oonposition  of  tbs  signal;  that  ia,  bow  auch  of  tba  energy  of  tbs  signal  is 
contained  in  each  different  frequency  band.  Finding  tbs  power  spectral 
density  can  be  thought  of  as  a  process  in  which  the  data  is  filtered  over  a 
very  narrow  frequency  band  (as  Af  ■*  0)  and  than  tbs  asan  square  value  of  that 
band  is  coaputed.  In  exact  fora: 

T 

Mf>  •  JF  t  |  x’ct.f.iftdt] 

O 

Often,  tbs  power  spectral  daasity  function  is  defined  as  tbs  Fourier 
trane fore  of  the  autocorrelation  function: 

+• 

G,(f)  "  2  l  R,(T)s’J2TTfTdT  , 

where  S* (t)  is  the  autocorrelation  defined  in  2.4. 1.3.  The  function  G«(f)  is 
a  poaitive,  real-valued  function. 

To  aotiaate  the  power  apectral  density  function,  the  date  were  analysed 
using  a  forward  fast  Fourier  tram  fora.  The  first  and  last  10%  of  each  data 
record  was  cos ins- tapered  to  avoid  the  Gibbs  phon canon  ( Benda t  and  Piersol 
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[1971]).  This  phenoBs non  ia  seen  aa  a  large  aplke  near  the  upper  and  lower 
cut-off  frequenciaa.  Thaae  large  overahoota  are  coaaon  near  digital  filter 
cut-offs,  and  the  Fourier  tranafora  behavea  aa  a  filter.  The  power  a pec t run 
eetinate  wea  co^uted  according  to: 

'vf  l*kl2  • 

where  Xk  ia  the  Fourier  coefficient,  h  ia  the  aaapling  interval,  and  N  is  the 
total  nunber  of  data  points. 


2. 4. 1.3  Autocorrelation  Function 

The  autocorrelation  function  describes  the  dependence  of  data  at  a  given 
point  ut  one  instant  in  tine  to  data  at  a  different  tine.  The  exact  autocor¬ 
relation  function,  for  a  given  neperation  tine  between  data  points,  is  given 
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x(t)  x(t+T)dt  . 
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The  quantity  R»(t)  is  real-valued  and  Bay  be  positive  or  negative.  R»(T)  has 
a  aaxiauB  at  nfl  and  is  syHotric  about  that  point,  i.a.,  R(t„)  =  R(-t0). 

The  autocorrelation  estiaate  is  coaputed  by  deteraining  the  power  spectral 
denaity  function,  and  then  coaputing  the  inverse  Fourier  tranafora  to  recover 
the  tiae  does in  and  the  autocorrelation  function.  This  eotinate  of  the 
autocorrelation  actually  yields  a  "circular”  function,  giving  spurious  results 
•t  larger  values  of  t.  However,  the  autocorrelations  of  the  well-pressure  and 
the  aass-flow  signals  die  out  very  quickly,  and  the  circular  part  of  the 
function  is  therefore  not  in  the  region  of  interest.  Finally,  the  auto- 


correlation  function  was  normalized  by  the  value  at  t  =  0,  giving  Rx(0)  =  1.0. 

2.4. 1.4  Croaa  Correlation  Function 

The  cross-correlation  function  describes  the  dependence  which  one  signal 
has  on  another,  as  a  function  of  the  delay  time  between  the  two  signals.  It  is 
determined  by  taking  the  average  product  of  x(t)  at  a  time  to  and  of  y(t)  at  a 
time  to  +  T  over  a  total  period  T,  as  T  approaches  infinity.  In  exact  fora: 

T 

Rxy(T>  j  X(t>  y(t+r)  dt  * 

0 

The  function  R*r(T)  is  a  real-valued  function  which  nay  be  positive  or 
negstive.  Unlike  the  autocorrelation  function,  it  does  not  necessarily  have  a 
maximum  at  t  =  0  nor  is  it  an  even  function.  However,  Rxy(T)  is  sjanwtric 
about  the  ordinate  when  x  and  y  are  interchanged;  i.e. ,  R*y (t)  =  Ry*(-T).  The 
value  of  the  cross-correlation  was  normalized  by  the  product  of  the  standard 
deviations  of  the  two  signals.  Thus,  if  x(t).and  y(t)  were  completely  similar 
at  a  specified  t,  then  R«y(T)  =  1;  likewise  if  H*r(T)  =  0,  then  the  two 
signals  are  said  to  be  independent  at  that  tine  delay. 

One  of  the  moat  important  applications  for  the  cross-correlation  function 
is  to  determine  the  time  required  for  a  signal  source  (that  is,  an  organized 
etructure)  to  pass  from  one  point  in  space  to  another.  If,  for  example,  two 
aignals  are  measured  in  a  flow  at  different  locations  and  a  peak  occurs  in  the 
crosscorrelation  function  at  a  certain  ?0(  it  can  be  said  that,  averaged  over 
time,  this  "event"  is  likely  to  take  T0  seconds  to  go  froa  one  point  to  the 
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other.  If  these  two  points  are  aligned  with  the  flow  direction,  a  convection 
velocity  for  this  event  can  be  found  by  taking  the  ratio  of  the  distance 
between  the  points  and  the  tine  delay  at  which  the  maximum  of  the  cross- 
correlation  occurs. 

'  To  estimate  the  cross-correlation  function,  two  different  Methods  were 
used:  fast  Fourier  transforms  and  a  direct  calculation  via  suns  and  products. 


Both  aethods  yielded  similar  results,  with  the  Fourier  nethod  having  values 
spproxinately  5*  higher,  due  in  part  to  a  circular  definition  of  the  cross- 
correlation  function.  Since  the  computation  tine  "for  the  Fourier  transform 
was  much  less,  it  was  used  almost  exclusively. 

As  discussed  previously,  a  10*  end  taper  was  applied  before  transforming 
the  data.  One  signal  was  packed  into  the  real  part  of  a  complex  aray,  the 


other  signal  into  the  imaginary  part,  and  then  a  forward  fast  Fourier  transform 
was  performed.  A  raw  cross-spectral  density  estimate  was  determined  as 
follows: 


vv  -  f  K  \i  • 


where  Yu  is  the  transform  of  y(t)  and  Xu  is  the  complex  conjugate  of  the 
transform  of  x(t).  Similarly  to  the  calculation  of  the  autocorrelation,  an 
inverse  Fourier  transform  of  Gxr(fk)  was  computed  and  normalized  by  the 


product  of  the  signal's  standard  deviations  to  determine  the  normalized 
cross-correlation  function,  R>r(t). 
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2.4. 1.5  Cross  Spectral  Density  Function 

Just  as  the  power  spectral  density  function  evolves  from  the  autocor¬ 
relation,  so  the  cross-spectral  density  function  evolves  froa  the  cross¬ 
correlation.  The  aost  interesting  aspect  of  the  cross-spectral  density  is  not 
the  cross-spectrum  but  rather  the  two  functions  which  result  froa  it:  the 
coherence  function  and  the  phase  function.  Since  the  cross-correlation  is  not 
an  even  function,  its  Fourier  transfers  (the  cross-spectral  density)  is  a 
coaplex  number: 


G„  (f)  -  C  (f)  -  J  Q  (f)  . 
xy  xy  J  xy 


The  real  part  of  the  cross-spectral  density  (the  coincident  spectral 
density  function)  can  be  thought  of  as  the  average  product  of  the  two  signals 
within  a  narrow  frequency  band,  divided  by  the  frequency  interval.  The 
imaginary  part  (the  quadrature  spectral  density  function)  is  the  sane  except 
that  either  input  signal  (but  not  both)  is  shifted  in  tiae  to  produce  a  90° 
phase  shift  at  the  frequency,  f.  In  polar  notation,  the  cross-spectral 
density  becomes: 


Vf)  ■  iVf),‘ 
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exy(f)  -  t«n_1[Qxy(f)/Cxy(f)] 
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The  coherence  function  is  then  defined  as: 


(f) 
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Physically,  the  coherence  function  tells  us  how  correlated  two  signals 
are  as  a  function  of  frequency.  This  is  an  extremely  valuable  tool,  since  it 
contains  correlation  values  for  each  different  frequency;  at  a  glance,  the 
frequencies  which  contribute  to  a  correlation  and  the  frequencies  which  die 
out  before  being  detected  at  the  second  location  can  be  determined. 

The  phase  angle,  9xy(f),  consitutes  the  phase  shift  between  the  two 
signals  at  a  particular  frequency.  From  the  phase  angle,  the  time  delay  in 
which  an  organized  event  moves  from  one  location  to  another  can  be  determined 
as  a  function  of  frequency: 


t (f )  -  6  (f)/2irf  . 

xy 

Furthermore,  if  the  separation  distance  between  measuring  locations  is 
known,  the  convection  velocity  can  be  determined  as  a  function  of  frequency 
also. 


2.4.2  Conditional  Sampling  and  Averaging  Techniques 

Conditional  sampling  and  averaging  was  utilized  to  better  determine  the 
existence  of  organized  structures  in  the  boundary  layer.  Conditional  sampling 
is  a  technique  under  which  special  consideration  is  given  to  segments  of  the 
signal  where  certain  conditions  are  fulfilled.  The  total  number  of  consecutive 
points  satisfying  these  conditions  make  up  an  "event".  It  is  believed  that 


these  "event*"  ere  associated  with  particular  fluid  dynamic  phenomena  and 
are  the  signature  of  characteristic  structures.  The  conditional  average  is 
then  defined  as  an  average  taken  over  all  the  events. 

There  are  two  types  of  conditional  sampling,  one-point  and  two-point.  In 
one-point  conditional  sampling,  the  same  signal  is  processed  twice;  once  to 
define  the  conditions  under  which  averaging  will  occur  (detection)  and  once  to 
actually  average  over  the  detected  events  (averaging).  In  two-point  condi¬ 
tional  sampling,  the  detection  and  averaging  are  done  upon  signals  acquired  at 
different  spatial  locations.  Two  point  conditional  sampling,  therefore,  allows 
time  and  length  scales  of  the  events  to  be  explored. 

For  this  experiment,  the  following  conditional  averages  were  done:  1) 
one-point  on  the  wall-pressure  signal,  2)  one-point  on  the  mass  flow  signal  at 
various  points  through  the  boundary  layer,  3)  two-point  on  the  mass  flow  (one 
wire  directly  above  the  other),  4)  two-point  with  the  wall-pressure  signal 
conditioned  on  the  mass  flow  at  different  points  (both  spanwise  and  stream- 
wise  separations) . 

Two  methods  of  conditional  averaging  were  utilized:  a  zero-crossing 
technique  and  the  variable-interval  time-averaging  (VITA)  method.  The  zero¬ 
crossing  technique  was  a  simple  first-attempt  to  find  a  typical  event.  Very 
briefly,  a  positive  event  was  detected  when  N  negative  data  points  were 
followed  by  N  positive  data  points  (and  vice-versa  for  negative  events).  The 
number  of  points,  N,  was  dependent  upon  how  stringent  a  detection  criteria  was 
desired.  In  this  series  of  tests,  N  varied  from  two  to  ten.  The  events  were 
all  sumsed  about  the  first  positive  point  (in  the  case  of  the  positive  event) 
and  then  averaged  over  the  number  of  detections. 
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The  VITA  technique,  first  developed  by  Blackwelder  and  Kaplan  [1976], 
uses  the  intermittent  character  of  the  short-time  variance  of  the  detection 
signal  to  detect  common  patterns  in  the  boundary  layer.  It  is  assumed  that 
peaks  in  the  short-time  variance  signal  correspond  to  regions  of  higher 
turbulent  energy  and  therefore  indicate  the  existence  of  organized  structures. 
The  short-time  variance  of  a  fluctuating  signal,  Q(xi,t),  is  defined  by: 

t+*r  t+»sT  2 

varCxi.t.T)  -  ^  J  Q2(xlts)ds  -  f  Q(x1,s)ds] 
t-%T  t-*sT 


where  T  is  the  averaging  period  for  the  short  time  variance.  As  the  variance 
period,  T,  goes  to  infinity,  the  short-time  approaches  the  long-time  mean 
square  level.  As  Johansson  and  Alfredsson  [1982]  noted,  there  is  a  close 
relation  between  the  short-time  variance  period  and  the  time-scales  contri¬ 
buting  to  the  short-time  variance.  A  type  of  band-pass  filter  behavior  is 
observed,  with  the  center  frequency  nearly  equal  to  the  inverse  of  the  short- 
time  variance  period. 

The  VITA  technique  assumes  that  an  event  occurs  when  the  short-time 

variance  exceeds  kO2  where  k  is  the  chosen  threshold  level.  The  reference 

nns 

point  for  the  conditional  averaging  is  taken  as  the  midpoint  of  the  event. 
Once  the  reference  points  have  been  calculated,  the  averaging  can  be  done  as 
follows: 

1  M 
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where  M  is  the  number  of  events,  tj  is  the  reference  point  for  the  individual 
event,  and  ¥  is  the  relative  tiae  from  the  reference  (so  that  the  signals  can 
be  examined  upstream  and  downstream  of  the  center  of  the  event).  The  events 
were  classified  as  either  positive  or  negative  depending  upon  the  slope  of  the 
signal  at  the  Midpoint  of  the  event. 

In  tne  VITA  technique,  two  important  parameters  can  be  varied  to  optimize 
results:  the  period  of  the  short-tiae  variance  and  the  detection  threshold. 
Optimization  implies  that  the  average  event  that  results  froa  this  process  is 
representative  of  frequently  occurring  events,  and  not  a  rare  event  or  an 
atypical  one.  For  each  type  of  data,  both  parameters  were  adjusted  so  that  the 
detection  criterion  would  be  "switched  on"  for  all  the  aajor  events  and 
"switched  off”  for  regions  of  lesser  activity.  These  Judgements  were  rather 
subjective  but  the  aajor  conclusions  drawn  froa  this  conditional  sanpling  were 
not  critically  affected  by  the  choice  of  sampling  criteria. 


PBBUMUU8T  BSSULT8:  BOOHDABY  LAYKB  OVIBVIXN 


3.1  Mean  Flow  Survey  of  the  Boundary  Layer 

Prior  to  taking  turbulence  neaaur— ante,  several  naan-flow  surveys  were 
taken  in  the  floor  boundary  layer.  One  of  the  aajor  goals  of  this  part  of  the 
experiment  was  to  deteraine  the  extent  of  the  two-dimensional ity  of  the  flow. 

Previous  surveys  of  the  tunnel  floor  boundary  layer  (Settles  [1975], 
Taylor  [1984])  illustrated  that  the  static  pressure  was  constant  across  the 
boundary  layer  in  the  y-direction.  To  neasure  the  spanwise  and  longitudinal 
variations  in  static  pressure,  a  four  inch  spanwise  (s-direction)  region  and  a 
twenty  inch  streanwise  (x-direction)  region  were  surveyed.  The  static  probe 
was  calibrated  againat  the  wall  static  pressure,  with  a  resulting  calibration 
coefficient  of  0.91.  The  spenwise  survey  was  taken  at  x  *  14.125  inches,  y  = 
1.625  inches.  As  can  be  seen  in  Figure  16,  the  maximal  peak-to-peak  variation 
waa  less  than  5k.  The  streanwise  survey  on  the  centerline  at  y  =  2.0  inches 
revealed  a  peak-to-peak  variation  of  6k  (corresponding  to  1.4k  variations  in 
Mach  number)  and  a  slight  adverse  pressure  gradient  (see  Figure  17).  The 
pressure  gradient  was  caused  by  the  growth  of  the  displacement  thickness,  and 
it  was  considered  small  enough  to  be  neglected. 

Pitot  pressure  profiles  (y-direction)  were  taken  at  five  spanwise  and 
nine  streanwise  locations;  a  typical  profile  is  seen  in  Figure  18.  The  pitot 
and  static  pressure  data  were  used  in  the  Rayleigh-Pitot  formula  (1)  to  derive 
the  Mach  number  profile.  A  4k  stagnation  temperature  gradient  was  assumed 
through  the  boundary  layer  (To  i,  =  1.04  T0oo ) .  This  was  shown  to  be  a  good 


approximation  for  this  boundary  layer  by  Taylor  [1984].  The  stagnation 
temperature  profile  was  uaad  along  with  the  Mach  mabar  profile  to  yield  the 
atatic  temperature  distribution  (2),  which  in  turn  gave  the  sonic  velocity 
(3).  Finally,  the  product  of  the  local  sonic  velocity  and  the  Mach  niaber 
resulted  in  the  local  flow  velocity.  All  velocity  profiles  were  found  to  be 
similar,  and  a  typical  profile  is  shown  in  Figure  19. 

(1)  P./*0 

(2)  T  -  To  (l+(y-l)/2  M2) 

(3)  a  -  /yRT 

Using  the  Van  Driest  [1951]  transformation,  the  velocity  profiles  were 
mapped  onto  the  incompressible  plane,  and  an  extensive  logarithmic  region  was 
found,  as  can  be  seen  from  the  typical  profile  reproduced  in  Figure  20. 

To  investigate  the  two-dimensionality,  two  spanwise  pitot  pressure  surveys 
were  taken  at  points  inside  the  boundary  layer.  These  profiles  were  taken  at 
the  same  x  position,  x  =  18.5  inches,  at  two  different  y  positions,  y  =  0.4  and 
0.15  inches.  A  10k  variation  in  Pitot  pressure  (see  Figure  21).  This  three- 
dimensionality  is  most  probably  due  to  the  constriction  of  the  side  wall 
boundary  layers  through  the  converging-diverging  nozzle;  however,  for  the 
purpose  of  studying  the  boundary- layer  structure,  these  effects  nay  not  be  too 
significant  and  they  were  discounted  in  this  study. 


{[(Y+D/2]  Mj2  }  Y/(Y 
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Preston  tube  aeasureaents  were  taken  et  43  streaawiae  locations  and  the 
skin- friction  waa  deduced  according  to  the  eethod  of  Hopkina  and  Keener  [1966]. 
All  naan  flow  reaulta  are  tabulated  in  the  Appendix. 


3.2  fluctuating  Wall-Pressure  Meeeurseente 

The  r.a.s.  wal 1-preesure  fluctuation,  when  nonaali zed  by  the  nean  wall 
static  pressure  of  3.34  psi,  was  0.015.  The  r.a.s.  can  also  be  non- 
diaensionalized  by  the  dynaaic  pressure  or  the  wall  shear  stress;  this  results 

r 

in  values  of  0.0025  and  2.2,  respectively.  These  values  are  approxiaately 
one-half  of  Kiatler  and  Chen's  [1963]  values  at  Mach  3,  but  are  in  good 
agreeaent  with  those  of  Speaker  and  Ailaan  [1966],  Chyu  and  Hanly  [1968],  and 
Coe  [1969].  While  it  is  true  that  high-pass  filtering  at  250  Hz  reduces  the 
r.a.s.  intensity  by  alaost  10*,  a  significant  portion  of  that  is  noise  assoc¬ 
iated  with  tunnel  operation.  The  low  values  are  probably  caused  by  the  lack  of 
resolution  at  high  frequencies  due  to  the  finite  size  of  the  transducers. 

According  to  Shewe  [1983],  the  "ideal"  transducer  (one  which  would  aeasure  100* 

2ruT 

of  the  r.a.s.  pressure)  has  a  non-diaensional  diaaeter,  — -  ,  of  20.  The 

vvall 

transducers  used  in  this  experiaent  have  a  non-diaensional  diaaeter  of  about 
200.  According  to  Shewe,  our  transducers  should  give  an  r.a.s.  pressure  which 
is  60*  of  that  given  by  the  ideal  transducer,  and  this  result  seesw  borne  out 
by  the  present  experiaent. 

As  aent.;oned  previously,  a  resonance  was  found  in  the  Kulite  spectra  at 
480  KHz.  In  addition  to  that  spike,  however,  other  abnoraalities  were  found 
in  the  spectra  all  the  way  down  to  40  KHz.  The  energy  distribution,  seen  in 
Figure  22,  shows  the  effect  of  these  transducer  resonances.  It  is  clear  that 
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the  ipactrw  beyond  40  KHz  ia  inaccurate.  Bowtvar,  Boat  of  tha  aoerfy 
appears  to  bo  contained  at  lower  frequencies  ao  tha  affact  of  thaaa  raeonancaa 
on  tha  r.n.s.  laval  ia  probably  aaall. 

The  probability  density  function  of  the  wall  praaaure  fluctuatiooa 
(Figure  14)  baa  a  Gauaaian  ahape,  aa  expect ad. 

3.3  Fluctuating  Mass-Flow  Heaaurwaanta 

The  r.a.a.  aass-f low  fluctuation  level,  noraalizad  by  the  local  aaaa  flow, 
for  a  typical  boundary- layer  traverse  ia  ahown  in  Figure  23  along  with  reaults 
froa  previous  work.  The  profiles  show  good  agreaaant  and  ware  very  repeatable. 
A  typical  power  spectral  density  plot,  taken  at  y/5  *  0.15  (Figure  24),  shows  a 
spectrua  without  any  large  apikes  or  abnonalitiee  well  past  tha  region  of 
interest.  Normalized  probability  density  distributions  for  different  locations 
within  the  boundary  layer  are  ahown  in  Figure  15.  Those  appear  aa  expected; 
with  the  edge  of  the  boundary  layer  approaching,  the  negative  skewness  indi¬ 
cates  an  increased  internittency  resulting  fron  one-sided  nass-flow  spikes  (or 
a  nass-flow  deficit  froa  the  Bean)  and  the  positive  skewness  near  the  wall 
indicates  internittency  of  the  opposite  sign. 

As  in  the  case  of  nulti-pressure  transducer  runs,  an  important  factor  in 
aulti-wire  runs  is  the  relative  phase  shift  between  the  channels.  The  hot-wire 
aneacaaters  were  all  identical  OISA  55M10  systsns  and  the  wires  were  only 
used  if  they  had  aiailar  frequency  responses.  Bence,  the  relative  phase 
shift  through  the  frequency  range  of  interest  was  ainiaized. 
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Length  and  tin  seal*  infornation  nm  determined  using  two-point  corre- 
1  at  ions  and  spectrin  functions.  Those  two-point  measurements  reveal  a  substan¬ 
tial  Mount  of  infornation  concerning  turbulent  ecales.  Nevertheless,  the 
results  are  atill  tins  averaged  and  they  are  not  indicative  of  the  instan¬ 
taneous  properties  of  the  flow  field. 

4.1  Pressure- Pressure  Results 

As  discussed  previously,  the  four  wall-pressure  transducers  were  nounted 
in  a  rotatable  plug  with  fixed  aeparation  diatancee,  €  ,  and  nsasursnents  were 
taken  with  a  variety  of  yaw  angles  between  the  flow  direction  and  the  line  of 
transducers.  Space-tine  correlations,  R(Si,  £t,  t),  are  shown  in  Figure  25 
for  three  different  yew  angles.  Two  trends  should  be  noted:  1)  the  aaxina  of 
the  apace- tiae  correlations  decrease  with  increasing  separation  distance,  and 
2)  the  naxiaa  of  the  space-tine  correlations  decrease  as  the  alignaent  of  the 
transducers  changes  froa  streanwise  to  spaowise. 

Broad-band  convection  velocities  were  deduced  froa  the  naxiaa  of  the 
•pace-tine  correlations  when  £a  *  0  (streanwise  allgnnsnt) ,  and  they  are 
•hown  along  with  previous  results  in  Figure  26.  The  trend  of  higher  con¬ 
vection  velocity  for  larger  transducer  separation  nay  be  explained  by  the 
Proposed  hierarchical  structure  of  the  boundary  layer.  This  states  that  higher 
frequency  notions  (snaller  scales)  are  found  closer  to  the  wall,  and  lower 
frequency  notions  (larger  scales)  are  found  toward  the  edge  of  the  boundary 


layer,  with  an  appropriate  diatributioo  of  acalaa  in-between.  Thus,  the 
larger  atructurea  connect  downstream  faatar  than  the  aaall-ecale  structures. 
Nhan  the  transducers  are  eoved  farther  apart,  only  the  large-scale  structures 
■ova  fast  to  be  detected  by  both  transducers  and  therefore  a  higher 

convection  velocity  is  computed. 

The  specs  correlation  pattern,  R(5i,  €a,  0),  of  the  wall-pressure  fluc¬ 
tuations  can  be  napped  out  as  a  set  of  isocorrelation  contours  as  in  Figure 
27.  The  results  conpare  favorably  to  those  obtained  by  Tan  at  al.  [1985]  in  a 
siailar  boundary  layer  and  by  Bull  [1967]  in  a  subsonic  flow.  These  contours 
give  an  indication  of  the  spatial  scales  of  the  pressure-producing  structures. 
Free  a  nearly  isotropic  pressure  field  at  small  a pacings,  the  transverse  scale 
quickly  becomes  such  sore  significant  than  the  longitudinal  scale. 

Figure  28  shows  a  set  of  isocorrelstion  ooetours  of  the  space-tine 
correlation  wax  ins,  IBx(Ci,  Ct,  t  ).  These  contours  give  inforaation  on 
the  decay  rates  of  the  pressure  producing  structures,  and  it  nay  be  seen  that 
the  decay  is  nuch  faster  in  the  transverse  direction  than  in  the  longitudinal. 

After  the  cross -spec t run  between  pressure  signals  was  conputed,  several 
related  functions  were  determined.  The  coherence,  phase  angle,  tine  delay, 
and  convection  speed  functions  are  shown  in  Figures  29,  30,  31,  and  32  for 
sero  yaw  angle.  The  coherence  function,  Figure  29,  indicates  that  the  higher 
frequency  pressure  structures  do  not  retain  their  shape  as  well,  or  for  as 
long,  aa  the  lower  frequency  (larger-scale)  pressure  structures.  This  trend  is 
expected,  since  the  larger  structures  span  a  greater  distance  than  the  high 
frequency  structures  and  are  readily  seep led  by  both  transducers.  Note  that 
the  lack  of  resolution  (limited  nunber  of  spectral  points)  is  responsible  for 


the  "peaky"  nature  of  the  coherence  function.  The  phase  angle  function,  Figure 
30,  gives  similar  results.  It  is  worth  noting  that  the  phase  angle  see—  to 
increase  linearly  over  the  range  from  0  to  40  KHz.  The  ti—  delay,  Figure  31, 
and  the  convection  speed,  Figure  32,  again  show  that  the  higher  frequency 
structures  are  moving  more  slowly  than  their  low  frequency  counterparts. 

4.2  Mass  Flow  -  Mm—  Flow  Results 

Two  hot  wires,  one  directly  above  the  other  at  a  fixed  separation  dis¬ 
tance,  were  traverse d  through  the  boundary  layer.  The  calculated  space-time 
correlations  at  different  distances  from  the  wall  are  shown  in  Figure  33.  The 
peak  values  of  the  correlations  are  quite  high,  reaching  a  maximum  of  0.65  near 
the  middle  of  the  boundary  layer.  More  importantly,  the  delay  time  corres¬ 
ponding  to  the  peak  of  the  space-time  correlation ,  Tmax  .decrease  a  from  20 
microseconds  at  the  floor  to  nearly  zero  at  the  edge  of  the  boundary  layer. 

The  value  of  Tmax  along  with  the  wire  — paration  distance  A  and  the  local 
mean  velocity  can  be  used  to  define  an  angle  6  .  From  the  diagrams  in  Figure 
34,  we  see  that: 
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The  angle  9  may  be  called  a  "structure  angle,"  in  that  it  is  associated 
with  an  average  large-scale  motion.  The  results  from  three  different  traverse s 
can  be  seen  in  Figure  35  as  a  function  of  position  in  the  boundary  layer.  The 
angle  is  small  near  the  floor,  increases  quickly  to  about  45*,  and  it  remains 


constant  at  this  value  throughout  70%  at  the  boundary  layer.  Finally,  the 
angle  shows  a  rapid  increase  at  the  edge  of  the  boundary  layer.  Note  that  the 
distribution  of  the  structure  angle  seems  to  be  independent  of  the  two  dif¬ 
ferent  separation  distances  chosen. 

The  distribution  of  0  is  in  accordance  with  Head  and  Bandyopadhyay’s 
[1981]  observations  in  a  subsonic  boundary  layer  at  low  Reynolds  numbers.  They 
observed  that  the  hairpin  vortices  displayed  small  angles  near  the  floor  and 
45*  through  the  central  portion  of  the  boundary  layer,  followed  by  a  slight 
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increase  near  the  edge. 

The  coherence  function  between  two  wires  (separated  by  0.16)  is  shown  in 
Figure  36  for  various  positions  in  the  boundary  layer.  Due  to  the  relatively 
large  separation  between  the  wires,  information  about  the  higher  frequency 
motions  in  the  flow  (above  approximately  50  KHz)  was  not  obtained.  Therefore 
Figure  36  decribee  the  behavior  of  the  low  frequency  (large-scale)  structures. 
Near  the  floor,  these  structures  do  not  have  a  high  coherence,  as  expected  from 
the  low  space-time  correlation  values  seen  near  the  floor  in  Figure  33.  How¬ 
ever,  from  y/6  =  0.2  to  the  edge  of  the  boundary  layer,  the  coherence  over  the 
low  frequency  range  is  quite  high  and  remarkably  similar  in  shape,  indicating 
that  they  scale  upon  the  boundary-layer  thickness. 

4.3  Pressure  -  Mass  Flow  Results 

Cross-correlations  were  computed  between  pressure  fluctuations  at  the 
wall  and  mass-flow  fluctuations  measured  at  various  points  within  the  boundary 
layer.  Figure  37a  shows  the  correlations  with  the  hot  wire  located  0.45  6 
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downstream  of  the  pressure  transducer ,  while  Figure  37  b  shows  them  for  a 
streamwise  separation  of  0.912. 

The  first  point  of  interest  is  the  rather  low  level  of  correlation,  with 
a  maximum  peak  of  0.22.  This  level  of  correlation  was  observed  for  even  the 
smallest  transducer  separations  and  can  be  ascribed  to  the  differences  in  the 
frequency  content  of  the  pressure  and  mass-flow  signals.  This  reduces  the 
correlations,  as  can  easily  be  demonstrated  by  calculating  the  correlation 
between  two  sine  waves  of  differing  frequencies. 

Secondly,  although  the  separation  of  the  transducers  is  doubled  between 
Figures  37a  and  37b,  the  level  of  the  correlation  remains  nearly  the  same. 
Furthermore,  the  general  shape  of  the  correlation  is  retained  (the  small  peak 
at  a  negative  value  of  f ,  followed  by  a  sharp  rise  to  the  major  peak).  Hence, 
the  structures  appear  to  retain  their  shape  and  coherence  as  they  are  convected 
downstream. 

Figures  37c  and  37d  show  correlations  for  ■  panwise  separations  between 
the  transducers  of  0.116  and  0.456  respectively  (both  also  have  a  streamwise 
separation  of  0.096).  We  see  that  Figure  37c  shows  a  slight  decrease  in 
correlation  level  when  compared  to  the  previous  figures  which  had  no  spanwiae 
separation.  In  addition,  in  Figure  37d  we  observe  that,  with  an  increased 
spanwise  displacement,  the  correlation  has  broken  down  completely,  suggesting 
that  the  structures  have  a  very  limited  spanwise  extent. 

4.4  Summary 

These  two-point,  time-averaged  methods  have  revealed  the  existence  of 
organized  motions  in  the  boundary  layer.  They  appear  to  scale  on  outer- 


50 


layer  variables  and  convect  downstream  a  significant  distance.  At  the.  same 
time,  they  have  a  very  limited  extent  in  the  transverse  direction.  Further¬ 
more,  those  average  structures  show  traits  which  are  similar  to  motions  found 
in  subsonic  boundary  layers. 
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Chapter  5 

CONDITIONAL  SAMPLING  RESULTS  AND  DISCUSSION 

5.1  Introduction 

The  time  averaged  measurements  given  in  Chapter  4  showed  the  presence  of 
45°  structures  in  a  high  Reynolds  number,  supersonic  boundary  layer  similar  to 
the  structures  observed  in  subsonic  flows  at  low  Reynolds  number.  To  inves¬ 
tigate  the  detailed  instantaneous  structure  of  these  organized  motions, 
conditional  sampling  techniques  were  used.  Two  different  schemes  were  used, 
and  the  details  were  given  in  Chapter  2.  The  measurements  were  confined  to  the 
outer  layer  for  two  fundamental  reasons:  1)  we  were  not  able  to  access  the  y4 
values  at  which  the  bursting  process  occurs,  since  the  minimum  y4  value  was 
about  1000,  and  2)  the  time  scale  needed  to  identify  the  bursting  process  (lOv 
/u*)  was  too  small  to  resolve  with  our  current  instrumentation. 

Although  the  shape  of  these  outer-layer  structures  and  their  frequency  of 
occurrence  strongly  suggests  that  the  bursting  process  is  intimately  involved 
in  their  development,  we  have  not  been  afforded  the  opportunity  to  try  to  link 
the  two  events.  Therefore  the  words  "event",  "organized  structure”,  organized 
motion",  "large-scale  structure",  etc.  will  refer  to  the  outer-layer  motions 
which  we  have  identified;  while  the  word  "burst"  will  refer  exclusively  to  the 
bursting  process  in  the  inner  portion  of  the  boundary  layer. 

6.2  One-Point  Pressure 

In  this  section,  a  single  record  length  of  pressure  data,  consisting  of 


25,000  points  (0.025  seconds),  was  examined.  The  results  are  virtually 
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independent  of  the  length  of  the  record  and  the  particular  record  chosen. 

The  wall-pressure  signal  was  initially  conditioned  upon  itself  using  the 
VITA  technique.  Here,  detection  and  sampling  occur  at  the  same  spatial 
location.  In  Figure  38,  a  typical  pressure  signal  is  shown,  along  with  its 
short-time  variance  with  a  period,  T,  of  30  microseconds.  When  scaled  with 
outer-layer  variables,  this  period  corresponds  to  a  non-dimensionalized 
"decision  time",  T*  =  TU„/6,  of  0.6.  It  appears  that,  with  an  appropriate 
threshold  value  k,  the  short-time  variance  signal  functions  quite  well  as  a 
detection  scheme;  it  is  "switched  on"  for  the  more  energetic  and  sudden  events 
and  "switched  off"  for  quiescent  periods  and  more  gradual  events. 

A  number  of  individual  positive  events  are  shown  in  Figure  39.  These 
events  were  detected  with  T*  =  0.2  and  k  =  0.8,  and  it  will  later  be  seen  that 
these  values  optimized  the  average  event.  The  similarity  of  the  events  is 
remarkable.  Note  that  the  "downside"  of  the  major  peak  in  most  events  does  not 
return  to  the  level  at  the  beginning  of  the  cycle.  That  is,  the  positive 
pressure  event  is  slightly  one-sided.  Also,  the  major  rise  in  pressure  is 
always  preceded  by  a  decrease,  whether  it  is  sharp  or  gradual.  These  two 
observations  are  reinforced  in  Figure  40,  which  shows  a  number  of  positive 
events  superimposed.  The  significant  features  are:  1)  the  very  sharp  rise 
about  which  the  event  is  centered,  2)  the  elevated  pressure  level  after  the 
cycle  is  over,  and  3)  the  dip  in  pressure  prior  to  the  steep  rise. 

After  identifying  the  individual  events,  an  ensemble  average  was  computed. 
The  threshold  dependence  is  seen  in  Figure  41.  The  major  features  of  the 
positive  event  seem  to  be  essentially  independent  of  threshold  value;  the  steep 
rise,  and  the  size  of  all  three  averages  are  approximately  the  same.  The  major 


difference  is  found  in  the  secondary  peak  on  the  downside  of  the  cycle,  where 
the  pressure  falls  off  as  the  threshold  decreases  and  the  number  of  events 
increases.  This  trend  is  due  to  the  summation  of  random  signals,  and  as  the 
number  of  events  increase  some  "washout"  occurs.  The  frequency  of  these 
events,  non-dimensionalized  with  outer-layer  variables,  f*  =  ffi/U^,  can  be 
seen  in  Table  2.  We  see  that,  as  the  threshold  increases  from  0.7  to  1.0,  the 
number  of  events  falls  off  by  one-half. 

The  effect  of  changing  the  period  of  the  decision  time  for  the  short-time 
variance  is  seen  in  Figure  42.  Again,  the  most  significant  portion  of  the 
average  event,  the  steep  rise,  is  virtually  unaffected  by  changing  the  detec¬ 
tion  parameters.  The  secondary  peaks,  however,  increase  with  the  length  of  the 
decision  time.  For  the  short-time  variance  to  remain  above  the  threshold  level 
with  an  increasing  decision  time,  there  must  be  activity  over  a  longer  period 
of  time,  and  the  individual  events  are  more  apt  to  have  an  additional  peak.  It 
is  interesting  to  note  that  as  the  decision  time  goes  from  0.4  to  1.0,  the 
number  of  events  stays  roughly  constant.  We  will  see  later  that  this  trend  is 
the  opposite  to  that  observed  for  the  mass-flow  events. 

The  threshold  and  the  decision  time  were  chosen  to  give  the  most  energetic 
average  event  while  retaining  a  significant  number  of  realizations.  The 
average  positive  event.  Figure  43,  and  the  average  negative  event,  Figure  44, 
were  detected  with  k  =  0.8  and  T*  r  0.2  (10  microseconds).  The  character  of 
these  average  events  is,  of  course,  similar  to  the  individual  events  from 
which  they  evolve.  Therefore,  we  see  the  characteristic  strong  rise  about 
which  the  event  is  centered,  the  decrease  in  pressure  before  and  after  this 
"wave  front",  and  the  presence  of  the  secondary  peaks. 
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From  peak-to-peak  the  positive  event  is  2.2  times  the  pressure  r.m.s. 
level,  while  the  negative  event  is  2.5  times  the  r.m.s.  level.  These  average 
events  are  significantly  larger  than  those  in  Figures  41  and  42  where  the 
detection  parameters  were  not  optimized.  The  non-dimensionalized  frequency  of 
occurrence  for  the  positive  events  is  0.16  and  for  the  negative  events  it  is 
0.21.  The  period  of  each  of  the  average  events  is  approximately  40  micro¬ 
seconds,  and  when  this  is  combined  with  their  frequency  of  occurrence  and  the 
length  of  the  total  record,  it  appears  that  these  typical  events  (positive  and 
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negative  combined)  occupy  about  25%  of  the  total  signal  time.  It  seems  that  we 
are  looking  at  a  relatively  common  event.  Considering  the  high  level  of 
activity  associated  with  the  events,  we  can  state  that  they  contain  the 
majority  of  the  turbulent  energy  in  the  flow.  Presumably,  they  will  also  carry 
most  of  the  turbulent  shear  stress. 

For  comparison,  a  conditional  average  using  a  different  detection  scheme, 
the  zero-crossing  technique,  was  also  performed.  In  this  case,  there  was  a 
single  detection  parameter  N,  where  N  is  the  number  of  consecutive  negative 
data  points  followed  by  the  same  number  of  positive  data  points.  Figure  45 
shows  the  average  positive  zero  crossing  (N  =  5,  8,  10)  and  Figure  46  shows  the 
average  negative  zero  crossing  (N  =  10).  The  effect  of  varying  N  can  be  seen 
in  Figure  45.  As  the  criterion  becomes  less  stringent  (N  decreasing)  the 
number  of  events  increase  and  the  secondary  peaks  drop  off  more  quickly.  Once 
again,  this  demonstrates  the  wash-out  due  to  the  summation  of  random  signals. 

As  with  the  VITA  technique,  however,  varying  the  detection  parameter  does  not 
effect  the  basic  character  of  the  average  event. 


Most  importantly,  the  zero-crossing  and  VITA  methods  give  approximately 
the  same  shape  for  the  average  events.  Also,  the  frequency  of  occurrence  is 
very  similar  to  that  deduced  using  VITA:  f*  =  0.17  for  the  positive  events 
and  f*  =  0.16  for  the  negative  events  (as  compared  to  0.16  and  0.21  with  the 
VITA  technique).  The  standard  deviation  from  the  event  frequency  was  computed 
to  be  0.21  for  the  positive  events.  That  the  standard  deviation  is  greater 
than  the  mean  implies  that  this  is  not  a  Gaussian  event;  rather,  the  distri¬ 
bution  of  the  frequency  of  events  is  one-sided,  biased  towards  higher  fre¬ 
quencies. 

r* 

Overall,  it  appears  that  the  wall-pressure  events  deduced  by  conditional 
sampling  do  not  depend  strongly  on  the  sampling  method  or  criteria,  and  it 
also  appears  that  the  events  are  reasonably  common  and  contain  significant 
turbulence  activity. 

5.3  One-Point  Mass  Flow 

The  VITA  technique  was  used  to  condition  the  streamwise  mass-flow  signal 
at  several  points  throughout  the  boundary  layer.  Initially  only  one  hot  wire 
was  used,  so  that  detection  and  sampling  occurred  at  the  same  spatial  location. 

Before  examinining  the  results,  we  need  to  consider  the  effect  of  the 
sampling  criteria  on  the  nature  of  the  deduced  positive  and  negative  mass-flow 
events. 

In  Figure  47  a  mass-flow  time  trace  and  its  corresponding  short-time 
variance  are  shown  for  a  y/6  value  of  0.03  (y*  =  350,  or  y  :  0.8  mm).  The 
decision  time  is  T*  =  0.2  (10  microseconds)  and  the  marked  threshold  is  0.8. 

Nearly  all  of  the  positive  events  which  were  detected  show  a  feature 


which  was  also  seen  on  the  positive  pressure  events;  namely,  once  the  event  is 
over,  the  velocity  is  at  a  higher  level  than  before  the  cycle  began.  This 
observation  is  reinforced  by  Figure  48,  which  shows  a  number  of  individual 
events  superimposed.  The  similarity  of  the  events  along  the  steep  velocity 
rise  and  the  uniformity  in  width  of  the  peaks  is  particularly  striking. 

Three  interesting  trends  appear  in  Figure  49,  which  shows  the  frequency 
of  positive  events  and  total  events  (positive  plus  negative)  against  position 
in  the  boundary  layer.  First,  as  one  might  expect,  the  frequency  of  events 
decreases  as  the  edge  of  the  boundary  layer  is  approached.  Lower  frequency 
(larger)  motions  are  expected  in  the  outer  regions  of  the  boundary  layer,  and 
therefore  fewer  events  will  be  detected.  Also,  no  events  occur  in  the  free- 
stream,  so  a  drop-off  is  imperative. 

Second,  the  number  of  total  events  and  positive  events  falls  off  drasti¬ 
cally  near  the  wall.  This  is  caused  by  the  band-pass  filter  characteristics 
of  the  short-time  variance,  with  the  inverse  of  the  short-time  variance  period 
being  the  center  frequency.  With  a  hierarchical  structure  of  frequencies  in 
the  boundary  layer  (as  discussed  in  Chapter  4),  fewer  events  will  be  detected 
near  the  floor  since  the  frequencies  of  the  most  prevalent  structures  will  be 
too  high  to  contribute  to  the  short-time  variance. 

Third,  the  total  events  near  the  floor  consist  almost  exclusively  of 
positive  events,  while  further  from  the  wall  the  split  between  positive  and 
negative  events  becomes  more  equal.  Again  referring  to  the  hierarchical 
structure  of  the  boundary  layer,  it  would  appear  that  the  negative  events  are 
slower  events;  i.e.,  they  have  a  lower  frequency  content.  Since  there  are 
fewer  low  frequency  structures  near  the  wall,  the  number  of  negative  events 
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detected  there  will  be  lese  than  that  detected  elsewhere.  This  explanation  is 
given  further  credence  by  Figure  50,  which  ahowa  the  frequency  of  total 
events,  positive  events,  and  negative  events  versus  the  short-time  variance 
period  for  the  station  nearest  the  floor.  We  see  that  the  number  of  negative 
events  increases  as  the  short-time  variance  period  increases.  Referring  to 
the  band-pass  filter  nature  of  this  technique,  we  can  say  that  as  the  detection 
scheme  becomes  exclusively  a  lower  frequency  detection  scheme,  the  number  of 
negative  events  increases.  Hence,  to  detect  more  representative  negative 
events  near  the  floor  (that  is,  more  common)  a  longer  short-time  variance 
period  should  be  used,  such  as  50  microseconds  (T*  =  1.0). 

If  we  examine  the  distribution  of  positive  events  in  Figure  50,  it 
appears  that  a  period  of  10  microseconds  (T*  =  0.2)  would  most  nearly  represent 
the  true  character  of  a  positive  event  near  the  floor.  Thus,  positive  events 
have  a  higher  frequency  content  than  negative  events.  It  appears  that  when 
using  the  VITA  technique,  different  decision  times  should  be  used  to  detect 
positive  and  negative  events.  This  result  is  in  contrast  to  the  great  simi¬ 
larity  between  positive  and  negative  wall-pressure  events.  Furthermore,  the 
detection  criteria  for  the  mass-flow  events  may  have  to  be  varied  from  point 
to  point  within  the  boundary  layer. 

Figure  51  displays  the  threshold  effect  at  jr/S  -  0.33  with  a  decision 
time  of  T*  s  0.2.  The  gradual  decrease  in  events  as  the  threshold  increases 
is  as  eroected. 


Figures  52  and  53  illustrate  the  optimum  conditionally  averaged  positive 
•nd  negative  mass-flow  events  for  three  positions  in  the  boundary  layer.  For 
both  the  positive  and  negative  events  a  threshold  of  0.8  was  acceptable.  The 


decision  times  were  based  upon  Figure  SO  for  y/6  =  0.03  and  were  chosen  to  be 
T*  =  0.2  for  the  positive  event  and  T*  =  1.0  for  the  negative  event.  Ideally, 
a  plot  of  occurrence-frequency  versus  decision  time  should  be  consulted  for 
each  point  in  the  boundary  layer. 

Table  3  ahowB  the  non-dimensionalised  occurrence  frequencies  for  various 
points  in  the  boundary  layer.  In  subsonic  flow,  Willmarth  and  Sharma  [1984] 
computed  f*  =  0.13  at  y/6  =  0.01;  which  compares  favorably  to  that  computed 
here  at  y/6  s  0.03. 

Now  that  the  question  of  sampling  criteria  has  been  discussed,  consider 
the  results.  From  Figures  52  and  S3,  it  is  clear  that  the  typical  positive 
mass-flow  event  is  stronger  and  more  energetic  than  the  typical  negative  event 
since  its  peak-to-peak  value  is  much  greater.  It  can  also  be  seen  that  the 
"wave  front"  of  the  positive  event  is  steeper  than  for  the  negative  event, 
implying  that  the  positive  event  contains  higher  frequency  motions. 

The  shape  of  the  average  positive  event  at  the  lowest  position  in  the 
boundary  layer  can  be  interpreted  as  an  organised  structure  moving  down  in  the 
boundary  layer  and  bringing  faster-moving  fluid  with  it.  The  shape  of  the 
positive  event  at  the  edge  of  the  boundary  layer  implies  that  faster-moving 
fluid  still  moves  down,  but  the  dip  prior  to  the  steep  rise  suggests  that  a 
large-scale  structure  containing  slower  fluid  is  moving  up  through  the  boundary 
layer. 

S.4  Three-Point  Mass  Flow 

Two-point  conditional  sampling  was  performed  on  the  data  taken  with  the 
three  hot-wire  probe.  The  three  wires  were  located  at  y/6  =  0.37,  0.46,  and 
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0.57  and  iha  detection  signal  was  the  short-time  variance  of  either  the  middle 
wire  or  the  lowest  wire.  When  an  event  was  detected,  all  three  signals  were 
sampled  simultaneously. 

The  three  mass-flow  signals  and  the  short-time  variance  of  the  lowest 
wire  are  shown  in  Figure  54  for  a  decision  time  of  T*  =  0.2  and  a  marked 
threshold  of  0.8.  The  similarity  of  the  three  mass-flow  signals,  especially 
in  the  proximity  of  detected  events,  is  remarkable  and  gives  strong  indication 
of  an  organized  structure.  The  inclination  of  these  structures,  previously 
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discussed  in  the  context  of  the  cross-correlation,  is  also  evident  here  by  the 
shifting  of  the  peaks  from  position  to  position.  The  organized  structures  and 
their  angles  of  inclination  are  even  more  evident  in  Figure  55.  Here  indivi¬ 
dual  events  from  all  three  channels,  conditioned  upon  the  middle  signal,  are 
plotted.  For  each  series  of  detected  events,  the  time  delay  of  the  peak 
between  positions  is  slightly  different,  implying  that  the  angle  of  inclination 
and  the  convection  velocity  vary  from  structure  to  structure. 

Figure  56  depicts  the  average  positive  events  from  the  three  wires,  all 
conditioned  upon  the  middle  wire  (k  =  0.8,  T*  =  0.2).  As  expected,  the 
average  event  from  the  conditioning  signal  is  much  more  distinct  than  those 
events  which  are  two-point  conditioned  averages.  However,  the  existence  of  an 
average  organized  structure  spanning  the  distance  between  the  three  wires  is 
clearly  indicated. 

A  structure  angle  cat.  also  be  determined  from  the  average  positive 
events  by  using  the  time  delay  between  the  points  where  each  event  crosses 
zero  during  its  steep  rise.  The  angles  which  result  (45*  from  bottom  to 
middle,  63*  from  middle  to  top,  and  51*  from  bottom  to  top)  are  very  similar  to 
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those  deduced  from  the  cross-correlation.  This  is  not  surprising,  and  it 
illustrates  that  the  deduced  structures  are  indeed  typical  of  the  boundary 
layer  and  important  in  characterising  it. 

5.5  Pressure-Mass  Flow  Two- Point  Conditioning 

Two-point  conditional  sampling  was  performed  on  the  data  taken  with  a 
single  normal  hot  wire  and  three  wall-pressure  transducers.  At  all  times,  the 
short-time  variance  of  the  mass  flow  signal  was  used  as  the  detection  signal. 

The  height  and  stream  wise  extent  of  the  organised  structures  was  investi¬ 
gated  by  placing  the  hot  wire  directly  downstream  of  the  three  wall-pressure 
transducers  (aligned  in  the  streamwise  direction)  at  four  different  y/6  values 
in  the  boundary  layer.  For  a  hot  wire/ pressure  transducer  separation  of  0.096, 
Figure  57a  shows  the  average  positive  conditioned  wall-pressure  events  for 
each  of  the  y/6  values  of  the  hot  wire.  Likewise,  Figure  57b  shows  the 
average  pressure  events  for  a  streamwise  separation  of  0.276  and  Figure  57c 
depicts  the  average  pressure  events  for  a  separation  of  0.456.  All  of  these 
pressure  events  were  conditioned  upon  the  mass  flow  (k  =  0.8,  T*  =  0.2)  and 
the  average  positive  mass-flow  svent  is  shown  for  sach  y/6  value  in  Figure  57d. 

The  first  feature  to  notice  is  the  excellent  sverage  pressure  event  which 
is  extracted  by  conditioning  upon  the  mass-flow  signal;  this  is  best  illus¬ 
trated  in  Figure  57a  for  a  y/6  value  of  0.03.  This  average  event  is  almost 
as  large  as  when  the  pressure  was  conditioned  upon  itself,  although  slightly 
different  in  shape.  The  difference  in  shape  is  most  probably  because,  in  the 
latter  caae,  all  pressure  events  were  summed  about  a  point  on  the  steep  rise, 
accentuating  that  part  of  the  event  and  neglecting  the  downside  of  the  peak. 
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Looking  at  Figures  57a,  b,  and  c  independently,  it  can  be  seen  that  while 
the  structure  deteriorates  slightly  as  y/6  increases,  the  strong  correlation 
between  events  at  different  y/6  values  indicates  an  organised  structure  of 
significant  height  (at  least  0.35).  When  the  three  figures  are  taken  together, 
we  see  that  over  a  stream  wise  distance  of  0.455  the  structures  have  retained 
much  of  their  original  coherence.  This  is  further  supported  by  Figure  58, 
which  shows  the  average  pressure  and  mass-flow  events  for  a  streamwise  sepa¬ 
ration  of  1.45.  At  this  distance,  a  lot  of  background  noise  has  appeared, 
but  the  organized  structure  has  clearly  preserved  its  basic  identity. 

The  •  panwise  extent  of  the  organized  structures  was  explored  by  intro¬ 
ducing  a  spanwise  separation  between  the  hot  wire  and  the  wall-pressure 
transducer  (in  addition  to  the  hot  wire  being  0.09  5  downstream)  and  measuring 
the  mass  flow  at  four  values  of  y/6.  Figure  59  shows  the  average  positive 
events  for  a  spanwise  separation  of  0.115  and  Figure  60  shows  them  for  a 
separation  of  0.236.  It  is  readily  seen  that  the  organized  structure  which 
retained  its  shape  in  the  streamwise  direction  has  lost  its  correlation  within 
a  very  short  spanwise  distance. 

5.6  Summary 

These  conditional  sampling  and  averaging  methods  have  reinforced  our 
ideas  from  Chapter  4  about  the  average  large-scale  structures  in  the  boundary 
layer.  We  have  seen  that  they  fill  the  height  of  the  boundary  layer,  they 
retain  their  shape  as  they  convect  downstream,  and  that  they  are  of  limited 
spanwise  extent.  The  instantaneous  character  of  this  technique  showed  that  the 
individual  events  are  very  frequent,  similar  in  shape  and  size  to  each  other 


and  to  the  average  event,  and  that  poeitive  events  are  much  more  energetic  than 
negative  events. 
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PINAL  DISCUSSION  AND  CONCLUSIONS 

..  -  -/  f'y-  ,  ■'  ^  • 

This  experimental  investigation  has  extended  our  knowledge  of  organized 
structures  in  turbulent  boundary  layers  by  examining  their  behavior  in  a 
compressible  flow.  Large-scale  structures,  inclined  at  45*  to  the  wall,  were 
found  to  fill  the  height  of  the  boundary  layer.  As  they  were  convected 
downstream,  these  structures  retained  a  great  deal  of  their  shape  and  charac- 
ter,  and  preserved  their  identity  for  at  least  1.5^ .  ^However,  the  spanwise 
extent  of  the  structures  was  limited.^  /  i  ,  ;  ,  j 

The  individual  events,  whether  pressure  or  mass  flow,  were  similar  to 
each  other  and  to  the  average  event,  and  they  were  very  energetic.  In  general 
the  positive  mass-flow  events  had  a  higher  frequency  content  and  were  larger 
in  amplitude  than  the  negative  events. 

It  would  appear  that  the  effect  of  compressibility  on  the  large-scale 
motions  is  very  small.  Certainly,  there  are  more  similarities  than  differences 
between  the  events  found  in  this  compressible  flow  and  those  observed  in 
incompressible  boundary  layers.  These  similarities  are  clearly  shown  by 
comparing  the  structure  angles  and  the  conditionally-averaged  results. 

For  example,  the  deduced  structure  angles  from  this  investigation  are 
consistent  with  Head  and  Bandyopadhyay’s  [1981]  observations  in  incompressible 
flow.  The  eddy  angle  is  low  near  the  floor,  increases  quickly  to  45°  (where 
it  remains  throughout  70%  of  the  boundary  layer),  and  increases  again  near  the 
•dge  of  the  boundary  layer. 
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Furthermore,  Johansson  and  Alfredsson  [1982]  used  VITA  conditional 
sampling  in  fully  developed  turbulent  water-channel  flow.  Figure  61  shows  a 
number  of  individual  positive  streamwise  velocity  events  from  their  flow. 
Qualitatively,  this  figure  compares  quite  favorably  to  Figure  48,  which  are 
individual  positive  mass-flow  events  from  this  investigation.  Both  figures 
show  the  individual  events  possessing  a  uniform  steep  rise  and  having  a  higher 
level  at  the  end  of  the  event  than  the  beginning.  In  addition,  the  average 
amplitude  of  the  events  in  both  cases  seems  to  be  approximately  four  times  the 
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r.m.s.  level. 

The  frequency  of  the  events  in  this  flow  also  compare  favorably  to 
incompressible  flow.  For  instance,  Willmarth  and  Sharma  [1984].  Normalized 
with  outer  flow  variables,  the  frequency  of  events  is  approximately  .15  close 
to  the  wall  (y/6  ~  .01). 

Several  examples  of  average  positive  velocity  events  in  incompressible 
flows  were  examined  and  found  to  be  strikingly  similar  to  the  mass  flow  events 
from  this  investigation  (Fig.  52).  Thomas  and  Bull  [1983]  performed  the  most 
extensive  measurements,  from  y/6  ~  .05  to  .91,  as  seen  in  Figure  62.  The 
character  of  the  average  event  changes  (with  increasing  height  in  the  boundary 
layer)  in  a  manner  similar  to  that  of  the  compressible  flow  investigated  here. 

As  position  in  the  boundary  layer  increases,  the  drop  prior  to  the  steep  rise 
(or  burst)  becomes  more  significant  in  both  cases. 

Johansson  and  Alfredsson  [1982]  and  Blackwelder  and  Kaplan  [1976]  both 
present  VITA  average  positive  events  from  positions  nearer  the  wall,  and  the 
results  are  reproduced  in  Figure  63  and  64,  respectively.  Despite  the  dif¬ 
ference  in  position  (y*  =  15  in  Figures  63  snd  64),  the  incompressible  events 
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are  nearly  identical  to  those  events  from  the  center  of  the  compressible 
boundary  layer.  The  amplitudes  of  the  incompressible  and  compressible  events 
are  also  very  similar. 

In  summary,  the  organized  structures  in  compressible  flow  have  a  large 
number  of  similarities  to  their  counterparts  in  incompressible  flow.  The 
large-scale  structures  are  similar  in  both  flows  and  more  surprisingly,  the 
incompressible  small-scale  structures  have  some  characteristics  of  the  compres¬ 
sible  large-scale  structures.  While  a  much  more  intensive  study  is  called 
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for,  this  preliminary  investigation  indicates  that  the  effect  of  compres¬ 
sibility  on  organized  structures  in  a  turbulent  boundary  layer  is  small. 
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APPEND I X  A 


HE AN  FLOW  DATA  TABULATION 


NOTE:  All  pressure*  arc  in  N/m**2  and  all  temperature*  are  ir,  Helvir. 


SPANUISE  STATIC  PRESSURE  PROFILE 
X-  39.  88  cm 
V«  4.  13  cm. 


Average  Stagnation  Preeiure*  6886E-O6 
Average  Stagnation  Temperature"  297.  2 


Z  (cm  ) 

Static  Preteure 

Z(cm.  ) 

Static  Pmstu 

9  080 

. 2268C-03 

1.  812 

2313E-G5 

9  083 

2269E-09 

1  743 

.  2297E-0S 

9  019 

.  2267E+09 

1  674 

2261E-05 

4  447 

2266E-09 

1  603 

2277E-03 

4  874 

. 2269E-03 

1  930 

2280E->0? 

4  802 

. 2261E+09 

1  498 

228 1 E-G3 

4  733 

.  2267E+09 

1  389 

2283E-03 

4  641 

226 IE *09 

1  317 

2284E-0? 

4  987 

.  2264E-03 

1  244 

22828-03 

4  918 

2266E-09 

1  174 

22B3E-03 

4  446 

2299E-09 

1  101 

2283E-03 

4  376 

.  2264E+09 

1  029 

2284E-03 

4  304 

.  2264E-05 

9982 

22B6E-03 

4  234 

2239E-03 

8867 

.  2283E-03 

4  160 

226 IE -09 

8163 

2289E-03 

4  086 

2291E-09 

7491 

2294E-05 

4  014 

2299E-09 

6717 

23OOE-03 

3  946 

2239E-03 

6027 

2304E-0? 

3  873 

2294E-03 

9304 

2319E-09 

3  802 

2269E-09 

4626 

2320E-03 

3  728 

2264E-09 

3896 

2319E-03 

3  696 

2273E-09 

3173 

2310E-03 

3  987 

2274E*09 

2494 

2300E-0S 

3  914 

22B0E-09 

1720 

2287E-G3 

3  442 

2279E-09 

1019 

2276E-05 

3  372 

2274E-09 

3296E-01 

2269E-C3 

3  301 

2272E-09 

0  OOO 

2269E-03 

3  229 

2268C-09 

-  6347E-01 

2299E-03 

3  199 

. 2280E-09 

-  1332 

2267E-03 

3  088 

2302E-09 

-  2077 

2275E-03 

3  014 

2307E-09 

-  2803 

2292E-03 

2  949 

2312E-09 

-  3900 

2309E-03 

2  871 

2307E*09 

-  4197 

2312E-C  3 

2  802 

2296E-09 

-  4938 

.  2320E-03 

2  730 

2293E-03 

-  9643 

231 9E-03 

2  697 

2286E-0S 

-  6369 

2304E-C3 

2  986 

22B2E-03 

-  7081 

2294E-0S 

2  919 

2277E-09 

-  7804 

2277E-03 

2  449 

2271E-09 

-  8490 

2269E-03 

2  374 

2274E-03 

-  9184 

2293E-03 

2  306 

2273E-09 

-  9932 

2232E-03 

2  238 

2280E-03 

-1  067 

22936-03 

2  169 

2286E-09 

-1  137 

2293E-03 

2  092 

2294E-09 

-1  208 

2233E-05 

2  022 

2309E-09 

-1  279 

2260E-03 

1  990 

231 4E-03 

-1  331 

2261E-03 

1  880 

23l9E*09 

-1  422 

2269E-03 

Zicm.  ) 


Static  Pritiun 


-1  494 

.  2282E*03 

-1  967 

. 2284E*09 

-1. 638 

2283E*09 

-1.  710 

22828*09 

-1  781 

22768*09 

-1.  892 

.  228084-09 

-1  922 

.  228784-09 

-1  994 

22948 ♦OS 

-2  068 

2304E*09 

-2.  139 

23038 *09 

-2  209 

22928 ♦OS 

>2  284 

22788*09 

-2  392 

.  22968*03 

-2  421 

22348*09 

-2  499 

.  22228*09 

-2  968 

22126*09 

-2  641 

.  22098*09 

-2  710 

.  22188*09 

-2  781 

.  22248*09 

-2  849 

22318*09 

-2.  920 

22308*09 

-2  992 

22398*09 

-3  069 

.  2244E*09 

-3  136 

.  22918*09 

-3  208 

.  22606*09 

-3  278 

22736*09 

-3  346 

22808*09 

-3  417 

.  22848*09 

-3  489 

.  22908*09 

-3.  999 

.  22998*03 

-3  634 

.  23088*09 

-3  703 

23298*09 

-3  774 

.  23478*09 

-3  843 

.  23628*03 

-3  914 

.  23698*03 

-3  986 

.  23638*09 

-4  099 

2348E*09 

-4  128 

2337E*09 

-4  202 

23268*09 

-4  276 

2324E*09 

-4  346 

.  23298*09 

-4  417 

.  23268*09 

-4  489 

23298*09 

-4  999 

23298*09 

-4  631 

23268*09 

-4  702 

2327E*05 

-4  773 

23298*09 

-4  847 

.  2330E*09 

-4  919 

23308*09 

-4  991 

23318*09 

-9  064 

23328*09 

-9  139 

23308*09 

SPANWISE  PITOT  PRESSURE  PROFILE 
tm  4*  44  CP 
V»  1  02  C* 

Avtragt  Stagnation  Pftiii/fl*  68806*06 
Avorago  Stagnation  To*poraturo»  294  0 


Z  (cp  ) 

9  080 
5  077 
9  020 
4  490 
4  874 
4  BO  3 
4  730 
4  *94 
4  986 
4  919 
4  444 
4  372 
4  300 
4  224 
4  194 
4  089 
4  019 
3  449 
3  879 
3  802 
3  731 
3  *61 
3  989 
3  916 
3  449 
3  374 
3  302 
3  231 
3  161 
3  089 
3  019 
2  948 
2  J79 
2  803 
2  739 
2  *62 
2  989 
2  916 
2  449 
2.  374 
2  303 
2  232 
2  162 
2  091 
2  022 
1  992 
1  881 
1  808 
1  739 
1  *68 
1  999 
1  926 


Pitot  Prottoro 


Z(cp  > 


1 9846*06 

1  499 

1987E*0* 

1  389 

1984E+06 

1  316 

19736*06 

1  242 

1987E+06 

1  168 

198 IE ♦06 

1  100 

1 986E*06 

1  028 

19786*06 

.  4962 

1989E*06 

8849 

1 9686*06 

8160 

1976E*06 

7414 

1562E*06 

6696 

1560E*06 

9990 

1 5*9E*06 

9272 

1 5976*06 

.  4933 

196 IE >06 

3608 

1 5726*06 

.  3116 

19626*06 

.  2410 

1 991E*06 

1700 

1 9336*06 

. 9679E-01 

1991E+06 

2321E-01 

1 9446*06 

3311 £-02 

1949E*06 

0  000 

19976*06 

-  4709E-01 

1549E*06 

-  *401 E 

-01 

19496*06 

-  1394 

1 9946*06 

-  2093 

1 5496*06 

-  2803 

1 9376*06 

-  3939 

19626*06 

-  4290 

15416*06 

-  4989 

1549E*06 

-  9688 

1 9*  IE  ♦06 

-  *420 

1 9336*06 

-  7100 

159*6*06 

-  7818 

19446*06 

-  8943 

1 530E*06 

-  4227 

19146*06 

-  4433 

1931E+06 

-1  066 

14476*06 

-1  142 

1930E*06 

-1  214 

1447E+06 

-1  282 

14896*06 

-1  396 

14496*06 

-1  42-* 

14796*06 

-1  497 

14*06*06 

-1  969 

14976*06 

-1  641 

. 14966*06 

-1  713 

1431E*06 

-1  789 

1 4306*06 

-1  894 

14386*06 

-1  929 

14496*06 

-1  498 

Pitut  Pmiute 

144iE*Cfc 
1427E*06 
1397E*0fc 
143iE*0o 
. 1435E*06 
. 1432E+06 
1444E*06 
1440E*06 
1«43E*06 
1439E*06 
14416*06 
14466*06 
1 453E*0& 
1493E*06 
1450E*06 
1494E-06 
. 1495E*06 
.  146BE+06 
14546*06 
1494E*06 
14*BE*06 
1 4 726*06 
1477E-06 
1477E*06 
1497E*06 
14B3E*06 
190lE*0o 
1 493E*06 
14B6E*06 
1506E*0fc 
1497E*06 
1926E*06 
19196*06 
1910E*06 
191 9E*06 
I 997E*06 
1939E*06 
1 949E*06 
1 959E*06 
1992E*06 
1944E*06 
1 939E*Oo 
1999E*06 
195lE*06 
1 5996*06 
1999E-06 
1 961E+06 
19716*06 
197*E*0fc 
1997E*0fc 
1564E*Ge 
196G£*0o 


Z  <  CM  > 


Pitot  Prtnuri 


-2  o** 

1 9796*06 

-2  142 

19648*06 

•2  214 

19728*06 

-2  283 

19648*06 

-2.  398 

.  19948*06 

-2  427 

19678*06 

-2  444 

19908*06 

-2  973 

19648*06 

-2  444 

19948*06 

-2  714 

.  19988*06 

-2  784 

19648*06 

-2  897 

19478*06 

-2  423 

19998*06 

-2  443 

19416*06 

-3  064 

19998*06 

-3  138 

19926*06 

-3  212 

19498*06 

-3  281 

19416*06 

-3  331 

19398*06 

-3  421 

19228*06 

-3  442 

19406*06 

-3  961 

19328*06 

-3  674 

19246*06 

-3  702 

19428*06 

-3  776 

19338*06 

-3  844 

19166*06 

-3  417 

19328*06 

-3  487 

.  19188*06 

-4  096 

.  19168*06 

-4  124 

19168*06 

-4  201 

19126*06 

-4  273 

19078*06 

-4  343 

19136*06 

•4  416 

19088*06 

•4  484 

19248*06 

-4  960 

13148*06 

-4  630 

19176*06 

-4  701 

19248*06 

-4  773 

19236*06 

-4  843 

19306*06 

-4  414 

1931E*06 

-4  484 

19608*06 

-9  062 

19698*06 

-3  133 

19878*06 

S'ANWISE  PITOT  PRESSURE  PROPILE 

X-  4*  99  eei 

v*  o  n  c* 


Average  Stagnation  Pretiure*  .  *8808*06 
Average  Stagnation  Temperature-  333  O 


Z<e«  > 

Pitot  Proaaure 

Z  t  cm  ) 

3  080 

1146E+06 

1.  432 

3  082 

11448+06 

1  383 

3  019 

1 1478*06 

1  311 

4  943 

11328*0* 

1  240 

4  874 

. 11398*0* 

1  171 

4  800 

. 11338*0* 

I  100 

4  738 

. 11498*06 

1  027 

4  *36 

11348*06 

9373 

4  383 

1 1 338*0* 

.  8843 

4  313 

11438*06 

8109 

4  441 

11398*06 

7410 

4  3*9 

114*8*06 

6*81 

4  299 

1 136E*06 

3993 

4  238 

11188*06 

327fc 

4  136 

11228*06 

4339 

4  089 

11178*06 

3833 

4  016 

11188*06 

3142 

3  943 

11218*06 

2429 

3  869 

111*8*06 

1723 

3  798 

11218*06 

9733E-01 

3  738 

11138*06 

3469E-01 

3  *37 

1 1 1 38*06 

0  OOO 

3  382 

11138*06 

-  63498-01 

3  313 

112*8*06 

-  13*9 

3  441 

11178*06 

-  2097 

3  369 

11248*0* 

-  2833 

3  398 

11398*0* 

-  3361 

3  234 

11238*0* 

-  4393 

3  133 

11208*0* 

-  4978 

3  083 

11228*06 

-  3736 

3  014 

11238*06 

-  *433 

2  940 

11348*06 

-  7133 

2  8*9 

11238*06 

-  7833 

2  797 

11318*06 

-  8361 

2  728 

11148*06 

-  9238 

2  *34 

11208*0* 

-  9939 

2  384 

11168*0* 

-l  069 

2  311 

11118*06 

-1  144 

3  439 

11018*06 

-1  217 

2  3*8 

10838*06 

-1  286 

2  393 

10828*06 

-1  338 

3  228 

10738*06 

-1  430 

3  138 

1082E*06 

-1  301 

2  090 

10698*06 

-1  374 

2  019 

106*8*06 

-  1  *43 

l  948 

10368*06 

-1  717 

1  876 

10498*06 

-1  788 

1  803 

lw 368 *06 

-1  839 

1  733 

10338*06 

-1  929 

1  *64 

1 043E*06 

-2  OOl 

1  393 

10308*06 

-2  072 

1  333 

10378*06 

-2  146 

Pitot  Pressuri- 

1041E*06 
1039E*06 
1033E*06 
.  1031E*0o 
10398*06 
1057E*06 
1037E*06 
104BE*06 
10S6E*Ofc 
1031E*06 
1061E*0i 
103SE*06 
103lE*0i 
1044E-06 
1056E*0c 
1047E*0o 
1044E*0i 
1030E*0o 
1039E-06 
1034E*0fc 
1033E*06 
I037£*0f 

1060€*0e 
1078E*0* 
1060E*0fc 
1038E*0e 
1064E-06 
I 066E*06 
1069E*0e 
10BSE*0c 
lOSTE-Ofe 
1093E*0b 
109*E*06 
lOPTE-Cfc 
I 101E*0o 
11028*06 
I 1038*06 
1099E*0e 
1 1088*0e 

ll07E*Ct 
1 1 148*0t 
1 1148*06 
1 1098*0e 
1 10oE*C>6 
1  1068*06 
1 1 1 4E*06 
1 l 068*06 
J i 1 6E-06 
1 !01E*06 
I  1  I08*0e 
1  1  0 JF*Oe 
1  106fc*06 


Z(c»  ) 

*ltot  Prtiturt 

-2  21* 

.  1 1 12E+06 

-2  284 

1101E+0* 

-2.  397 

.  1116E+0* 

-2  430 

10**E«-06 

-2.  900 

1100C4-0* 

-2  974 

1 107E+O* 

-2  649 

1101E+06 

-2  71* 

1 104E+06 

-2  789 

10*1E«06 

-2  89* 

10*8E+O* 

-2  *24 

10*9E*06 

-2  *97 

10*7E*0* 

-3  0*7 

10**£«-O* 

-3  138 

108*E*06 

-3  213 

10*1E*06 

-3  282 

1083E*0* 

-3. 391 

10*7E<-06 

-3  422 

10*9C»06 

-3  4*1 

10*9E»0* 

-3  9*4 

10*9E«06 

-3  *2* 

10*9€*0* 

-3  703 

. 1088E+O* 

-3  779 

1079E*0* 

-3  843 

10886*06 

-3  *14 

107BE*06 

-3  *86 

1083E*06 

-4  097 

.  ioao£«-o* 

-4  128 

107*E*O6 

-4  200 

10806*06 

-4  2** 

1078E*0* 

-4  343 

10746*06 

-4  414 

. 10*16*04 

-4  48* 

. 1041E*0* 

-4  998 

1072E*06 

-4  *28 

10**E*0* 

-4  *** 

1074E*0* 

-4  771 

10806*06 

-4  841 

1087E*0* 

-4  *1* 

. 10*86*04 

-4  *8* 

10*26*06 

-9  0*1 

. 11026*0* 

-9  134 

111 1E*0* 

VERTICAL  PITOT  PRESSURE  PROFILE  AT  A  SPANWISE  STATION 
»■  44  99  cm 
2-  3  02  cm 


Av*r«|« 

6874E*06 

Avaraga 

-  260  9 

Y<«m  ) 

Pit*>t  Pr*«»gri 

Static  Prttiun 

U(M/« ) 

2540E-01 

60976*05 

23066*05 

366  4 

7A81E-01 

73246*05 

2306E*05 

399  1 

1222 

8347E*0S 

2306E*05 

421  1 

1725 

92366*05 

23066*05 

437  9 

2212 

98366*05 

23066*05 

447  9 

2674 

10346*06 

2306E*0S 

455  6 

32ie 

10836*06 

23066*05 

462  6 

3672 

111 26*06 

23066*05 

466  8 

4136 

11616*06 

23066*05 

473  0 

4654 

1191E*06 

23066*05 

477  0 

5089 

12246*06 

23066*05 

480  9 

6045 

12926*06 

23066*05 

489  5 

7066 

13586 <-06 

23066*05 

496  7 

8026 

14206*06 

2306E*0S 

502  7 

8971 

14796*06 

23066*05 

507  7 

9901 

15336*06 

. 23066*05 

512  5 

1  090 

.  10916*06 

2306E*0S 

518  0 

1  189 

.  16736*06 

23066*05 

524  0 

1  285 

17386*06 

23066*05 

529  1 

1  381 

17986*06 

23066*05 

533  2 

1  476 

18636*06 

2306E*03 

537  2 

1  67] 

.  2OO4E*06 

2306E*05 

546  1 

1  867 

21086*06 

23066*05 

551  9 

2  059 

2251E*06 

23066*05 

559  2 

2  251 

23406*06 

23066*05 

562  5 

2  443 

24206*06 

. 23066*05 

565  7 

2  637 

24626*06 

23066*05 

567  2 

2  833 

24866*06 

23066*05 

568  3 

3  028 

24966*06 

23066*05 

568  8 

3  218 

.  24966*06 

. 23066*05 

568  8 

3  414 

.  2498E*06 

23066*05 

568  8 

3  613 

24966*06 

23066*05 

568  8 

3  807 

251 1E+06 

. 23066*05 

569  4 

4  001 

25026*06 

23066*05 

568  8 

4  191 

25136*06 

23066*05 

569  4 

4  389 

24996*06 

23066*05 

566  8 

4  583 

25306*06 

2306E*0S 

570  5 

4  780 

25206*06 

230  6*05 

569  9 

4  974 

25306*06 

23066*05 

570  5 

5  169 

25246*06 

23066*05 

569  9 

(1 


I  276 
I  437 
I  554 
I  650 
1  7  J  C 
X  759 
1  804 

1  832 
X  874 
X  903 
X  931 
X  995 

2  051 
2  101 
2  143 
2  186 
2  235 
2  292 
2  342 
2  364 
2  427 
2  526 
2  597 
2  689 
2  738 
2  788 
2  816 
2  830 
2  83-' 
2  837 
2  63' 
2  837 
2  844 
2  837 
2  844 
2  837 
2  859 
2  851 
2  859 
2  851 


m  *  *  *  *  »  u  u  u  u  u  u  w  pj  pj  r ;  pj 


I'Li'ti'lj 


Li 'I.iI.i  Ia»* Lt ' 4 


VERTICAL  PITOT  PRESSURE  PROFILE  AT  A  SPANUISE  STATION 
X*  46.  99  cm 
2 *  1.  91  cm. 


Average  Stagnation  Pressure*  .  68736*06 
Average  Stagnation  Temperature*  297  7 

V<cm  )  Pitot  Pressure  Static  Pressure  U(m/s>  n 


2940E-01 

76B1E-01 

1218 

1723 

2208 

2677 

3218 

3692 

4166 

4631 

9089 

6071 

7047 

0007 

8948 

9908 

1  091 
1  189 
1  284 
1  379 
1  476 
I  669 
1  367 
060 
293 
443 
632 
832 

029 
220 
419 
608 
809 
003 
196 
397 
980 
780 
976 
171 


.  9896E*09 
. 69996*09 
. 78906*09 
. 89696*09 
.  91066*09 
9917E*09 
.  99606+09 
. 10436*06 
10706*06 
. 10966*06 
1 124E+0A 
1197E+06 
. 12746*06 
13306*06 
1361E*06 
1492E+06 
.  1494E*06 
19B3E*06 
1618E*06 
1673E+06 
. 17486*06 
.  IB62E+06 
. 1969E*06 
.  2090E+06 
.  2297E+06 
23088*06 
. 2399E*06 
.  24346*06 
.  2469E *06 
2490E*06 
24896*06 
29046*06 
2916E*06 
2914E*06 
.  290BE*06 
291 8E*06 
. 2910E*06 
2939E*06 
2930E*06 
. 2999E+06 


2280E+09 

398.  8 

1.  293 

22B0E+09 

389  7 

1.  402 

2280E+09 

410  7 

1  912 

22806*09 

424.  7 

1  989 

22806*09 

434  9 

1  646 

22B0E*o3 

441  9 

1  689 

22806*09 

448.  9 

1.  739 

22806*09 

496  3 

1.  7B2 

22B0E*09 

499  4 

1.  804 

22806*09 

463  6 

1  832 

22806*09 

466.  6 

1  893 

2280E*O9 

476  6 

1  924 

22806*09 

489.  1 

1.  988 

22806*09 

491.  4 

2.  037 

2280E*09 

494  8 

2.  066 

2280E+09 

903.  4 

2.  136 

2280E*09 

907  3 

2  172 

2280E*09 

919.  2 

2.  243 

22806*09 

917  3 

2.  264 

22B0E*09 

921  6 

2  306 

22806*09 

927.  3 

2.  363 

2280E*09 

934  9 

2  441 

22806*09 

940  6 

2  912 

22806*09 

947  7 

2.  997 

22806*09 

996  1 

2.  703 

2280E*09 

998  2 

2  738 

22806*09 

961.  9 

2.  788 

22806*09 

963.  7 

2.  816 

22806*09 

969  3 

2.  837 

22806*09 

969.  9 

2.  844 

2280E+09 

969.  9 

2  844 

22806*09 

967.  0 

2.  899 

2280E*09 

967  9 

2.  866 

22806*09 

967.  9 

2.  866 

2280E*09 

967  0 

2.  899 

22806*09 

967  9 

2.  866 

2280E*09 

967  0 

2.  899 

22B0E+09 

968  0 

2.  873 

2280E+09 

968  0 

2.  873 

22806*09 

969.  1 

2.  887 
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VERTICAL  PITOT  PRESSURE  PROFILE  AT  A  SPANWISE  STATION 
X-  46  99  cm 
2*  0  0  c* 


Avfragi 

Stagnation  Prtuuri*  . 

6873E+06 

Avingt 

Stagnation  Tamparatura 

-  251. 0 

V <  cm  > 

Pitot  Praaaura 

Static  Praaaura 

U<m/*> 

H 

3334E-01 

. 5916E+05 

.  2252E+05 

358  6 

1.  274 

8739E-01 

. 71B6E+05 

.  2252E+05 

392  1 

1  441 

.  1297 

. 8046E+05 

. 2252E+05 

410  5 

1  540 

.  1819 

. B738E+05 

. 2252E+05 

424  1 

1  618 

2299 

. 9368E+05 

2232E+05 

435  2 

1  685 

2734 

. 9813E+0S 

22S2E+05 

442  6 

1  731 

3304 

. 1024E+06 

. 2252E+05 

449  4 

1  775 

.  3735 

1050E+06 

.  2252E+05 

452  5 

1  797 

.  4238 

.  1090E+-06 

.  2252E+05 

458  7 

1  839 

4711 

1 130E+06 

. 2252E+05 

463  S 

1  874 

9172 

1 160E+06 

. 2252E+05 

467  7 

1  903 

6166 

1225E+06 

2252E+05 

475  3 

1  95c 

7156 

.  1286E+06 

. 2252E+05 

482  5 

2  016 

8090 

1363E+06 

. 2252E+05 

490  4 

2  080 

9020 

1399E+06 

. 2252E+05 

493  6 

2.  108 

9991 

.  148lE-*-06 

. 2252E+05 

501  0 

2  172 

1  101 

1557E+06 

. 2252E+05 

508  0 

2.  235 

1  197 

.  1607E+06 

. 2232E+03 

51  1  6 

2  271 

1  292 

.  1675E+06 

. 2252E+05 

517  4 

2  32e 

1  386 

.  1741E+06 

.  2252E+05 

521  5 

2  370 

1  483 

.  1792E+06 

. 2252E+05 

525  5 

2  412 

1  676 

.  1915E+06 

.  2252E+05 

533  0 

2  497 

1  873 

.  2014E*06 

. 2252E+05 

538  1 

2  561 

2  070 

.  2163E+06 

. 2252E+05 

546  0 

2  660 

2  262 

.  2275E+06 

. 2252E+05 

551  1 

2  731 

2  450 

.  2342E+06 

. 2252E+05 

553  7 

2.  774 

2  641 

.  2374E+06 

. 2252E+05 

554  7 

2.  795 

2  839 

.  2439E+06 

. 2252E+05 

557  9 

2  837 

3  038 

. 2450E+06 

. 2252E+05 

558  5 

2  844 

3  329 

.  2477E+06 

. 2252E+05 

559.  5 

2  85® 

3  423 

. 2465E+06 

2252E+05 

559  0 

2.  851 

3  618 

. 2498E+06 

.  2252E+05 

560  6 

2  873 

3  812 

. 2491E+06 

. 2252E+05 

560  1 

2  866 

4  010 

. 2515E+06 

.  22S2E+03 

561  1 

2.  880 

4  202 

. 2505E+06 

. 2252E+05 

560  6 

2  e73 

4  393 

.  2516E  +  06 

2252E+05 

561  1 

2  880 

4  586 

. 2509E+06 

2252E+05 

561  1 

2  880 

4  787 

2523E+06 

2252E+05 

561  6 

2.  887 

4  985 

2513E+06 

2252E+05 

561  1 

2.  880 

5  178 

2523E+06 

2252E+05 

561  l. 

2  887 

VERTICAL  PITOT  PRESSURE  PROFILE  AT  A  SPANWISE  STATION 
X-  46.  99  cm 
Z-  -1.  31  cm. 


Average  Stagnation  Preeeure-  6S77Er06 
Average  Stagnation  Temperature-  237.  0 


Y<  cm.  ) 

Static  Pressure 

U(m/ft ) 

H 

2540E-01 

.  6396E*03 

.  2281E+03 

374.  1 

1  327 

7000E-01 

.  7462E+03 

.  2281 ErOS 

401.  0 

1.  462 

1191 

842SEr03 

. 2281E+05 

421  2 

1  S72 

1679 

.  9172E+05 

. 2281E+0S 

435.  3 

1.  633 

2123 

981 3ErOS 

. 2281 Er03 

443  7 

1.  717 

2669 

1038E+06 

.  228 IE +03 

4S4  9 

1  779 

.  3119 

1074E+O6 

. 2281 ErOS 

460  3 

1  811 

.  3641 

1 120E+06 

. 2281 ErOS 

466  S 

1  833 

4117 

1 1 31E+06 

. 2281E+03 

470.  S 

1.  881 

4337 

1188Er06 

. 2281 ErOS 

475.  S 

1  917 

5077 

.  1237E+06 

. 2281Er03 

481.  3 

1  959 

.  6022 

1297E+06 

.  2281 ErOS 

487.  7 

2.  009 

o937 

1 366E+06 

.  2281 ErOS 

494.  9 

2  066 

7920 

.  1407E+06 

. 2281 ErOS 

499.  0 

2.  101 

3944 

.  1490E+06 

. 2281 E+05 

S06.  5 

2.  16S 

990S 

.  1363E+06 

2281 ErOS 

513  7 

2.  228 

1.  OSS 

1629E*06 

.  2281 ErOS 

SIB  1 

2.  271 

1  179 

1709E+06 

.  2281 ErOS 

524  8 

2  335 

1.  276 

1771E*06 

.  2281 E+OS 

528  9 

2.  377 

1  378 

1824E+06 

. 2281 ErOS 

332.  9 

2  420 

1.  475 

.  1933E-KJ6 

. 2281 ErOS 

339  S 

2.  490 

1  667 

.  2038E*06 

.  2281 ErOS 

S4S.  4 

2.  361 

1  838 

. 2173E+06 

.  2281 ErOS 

332  8 

2  693 

2  051 

2277E+06 

.  2281 ErOS 

557  * 

2  717 

2.  249 

. 241 SE+06 

2281 ErOS 

363  4 

2.  802 

2.  443 

2453E+06 

.  2281 ErOS 

364  2 

2  823 

2  632 

2302E+06 

.  2281 E+OS 

563.  7 

2.  831 

2.  824 

. 2493E+06 

.  2281 ErOS 

565  7 

2  851 

3  019 

2912E+06 

.  2281 ErOS 

366  2 

2  839 

3  213 

2499E+06 

.  2281 ErOS 

365.  7 

2.  BS1 

3.  414 

. 2513E+06 

.  2281Er05 

366  2 

2.  859 

3  608 

. 2302E+06 

2281 ErOS 

56  S'.  7 

2  851 

3.  796 

. 2310E+06 

.  2281 ErOS 

366  2 

2.  859 

3  993 

2498E+06 

.  2281E+05 

36S.  7 

2  831 

4  193 

.  2314E+06 

.  2281 Eros 

366  7 

2  866 

4  387 

. 290SE+06 

. 2281 Eros 

366  2 

2  839 

4  373 

. 2524E+06 

228 1 ErOS 

566  7 

2.  866 

4  775 

. 2329E+06 

.  22B1E+0S 

367  2 

2  873 

4  972 

2348E+06 

.  2281 ErOS 

367  8 

2  880 

3  166 

2349E+06 

2281 ErOS 

567  8 

2  880 

VERTICAL  PITOT  PRESSURE  PROFILE  AT  A  SPANWISE  STATION 
X»  46  99  e* 

Z-  -3.  04  cm 

Avaraga  Stagnation  Praaaura*  .  6B77E+06 
Avaraga  Stagnation  Tamparatura"  250.  9 

Y ( cm. )  Pitot  Praatuna  Static  Praaaura  U(m/*> 


2351E-01 

. 6054E+05 

.  7945E-01 

. 7436E+05 

1195 

. 8442E+05 

.  1717 

9195E+05 

2197 

. 9750E+05 

.  2651 

. 1017E+06 

3199 

.  106OE»06 

.  3622 

.  1 1 16E+G6 

.  4159 

.  1147E+06 

.  4609 

.  117S£*06 

.  5058 

.  1212E+06 

.  6060 

1293E+06 

7035 

1354E+06 

.  7973 

. 1417E+06 

.  8922 

.  1484E-*06 

9B09 

.  1546E+06 

1  009 

.  1616E+06 

1.  107 

. 1691E+06 

1  281 

. 1762E+06 

1  377 

.  1B41E+06 

1  473 

.  1899E-»06 

1  666 

. 2045E+06 

1  064 

. 2175E+06 

2  058 

. 2344E+06 

2  252 

. 2430E+06 

2  440 

.  2473E+06 

2.  632 

.  2503E+06 

2  829 

.  2500E+06 

3  028 

.  2506E*06 

3.  217 

.  2496E+06 

3  414 

.  2484E+06 

3  607 

.  249BE+06 

3.  804 

.  2502E+06 

4  001 

.  2522E+06 

4  192 

.  2530E+06 

4  384 

.  2525E+06 

4  580 

2544E+06 

4  778 

2553E*06 

4  975 

2567E+06 

5  167 

.  2567E+06 

2239E+05 

363  1 

1  296 

2239E«-05 

398  8 

1  476 

2239E+05 

419  9 

1  593 

2239E+05 

433  6 

1  674 

2239E+05 

442  7 

1  731 

2239E+05 

449  5 

1  775 

2239E+05 

455  9 

1  818 

2239E+05 

463.  0 

1  8o7 

2239E+05 

466  9 

1.  896 

2239E+05 

470  B 

1  924 

2239E+05 

474  6 

1  952 

2239E*05 

483  8 

2  023 

2239E+05 

490  7 

2  060 

2239E+05 

496  5 

2.  129 

2239E+05 

503  0 

2  186 

2239E+05 

508  3 

2.  235 

2239E+OS 

513  5 

2  285 

2239E+05 

519  2 

2  342 

2239E*05 

524  7 

2.  393 

2239E+05 

530  0 

2  455 

2239E+05 

533  0 

2  490 

2239E+05 

541  4 

2.  589 

2239E+05 

548  6 

2.  682 

2239E+05 

556.  3 

2  78S 

2239E+05 

559  3 

2.  837 

2239E+05 

560  6 

2  866 

2239E+05 

561 .5 

2  887 

2239E*05 

561  0 

2.  880 

2239E+05 

561  5 

2.  887 

2239E+05 

561.  0 

2.  080 

2239E+05 

560  5 

2  87? 

2239E«OS 

561  0 

2.  8SC 

2239E+05 

561  0 

2.  8S0 

2239E+05 

562  0 

2.  894 

2239E+05 

562  5 

2  901 

2239E+05 

562  0 

2  894 

2239E+05 

563  0 

2  908 

2239E+05 

563  6 

2  915 

2239E+05 

564  1 

2.  922 

2239E+05 

564  1 

2.  922 

STREAMWISE  STATIC  PRESSURE  PROFILE 
Y*  3  08  cm 
2«  0  0  cm 


Average  Stagnation  Preiture*  6B90E+06 
Average  Stagnation  Temperature^  234.0 


X  ( cm.  ) 

Static  Prtiture 

X  <  cm.  > 

Static  Pre*< 

16.  31 

.  2198E+03 

19  90 

. 2214E+03 

16.  S3 

.  220 IE *03 

19  96 

2214E+03 

16  59 

.  2205E+03 

20.  02 

. 2214E+05 

16  73 

2199E*03 

20.  09 

.  2222E+03 

16.  81 

.  2196E+03 

20.  16 

. 2220E+03 

16.  86 

.  219BE+03 

20  22 

. 2222E*03 

16.  92 

.  2202E*03 

20  29 

. 2223E+03 

17  00 

.  2206E+03 

20.  33 

.  2217E+03 

17  06 

.  2203E+05 

20.  42 

.  2214E+03 

17.  13 

.  2204E+03 

20  49 

.  221 BE *03 

17  20 

.  2201E*03 

20  33 

.  221 1E*03 

17.  26 

.  219SE+03 

20.  62 

.  2212E*03 

17.  32 

2199E+05 

20.  69 

.  2216E*03 

17  40 

2200E*03 

20  76 

.  221 3E*05 

17  46 

.  2202E+05 

20  81 

2217E*03 

17.  32 

2206E+05 

20  90 

.  2212E+03 

17.  39 

.  2204E*03 

20.  93 

. 221 1E*0S 

17  66 

.  2202E+03 

21.  01 

. 221 1E*03 

17  72 

.  2201E+03 

21  08 

.  2212E+03 

17.  78 

2197E+03 

21  13 

22l0E*03 

17  83 

.  2196E*05 

21  22 

.  2209E*03 

17  92 

.  2199E+03 

21  28 

.  22 1 3E+03 

17.  96 

.  2200E*05 

21  33 

. 2213E*03 

18.  04 

. 2200E+03 

21  40 

. 2208E+03 

18  11 

. 2206E+05 

21  47 

.  22 1 0E*03 

18  17 

.  220 IE *03 

21  33 

. 2210E+03 

18  25 

2202E*03 

21  61 

. 2203E*03 

18  30 

.  2202E+05 

21  67 

. 2212E+03 

13  37 

.  2201E*05 

21.  74 

. 2202E+05 

13  44 

.  2202E+0S 

21.  81 

2207E+03 

18.  51 

.  2205E+05 

21.  88 

.  221 1E*03 

18  37 

.  2204E+03 

21  94 

. 2206E*03 

13  63 

2199E+03 

22.  02 

. 2214E+03 

18  71 

.  2197E*03 

22.  09 

.  2219E*03 

18  76 

. 2193E+03 

22.  14 

.  2226E*03 

18  84 

.  2199E*03 

22.  21 

.  2227E*03 

18  91 

.  2207E*03 

22.  28 

2230E+03 

18  93 

. 2210E+03 

22  33 

. 2224E+03 

19.  03 

.  2218E*0S 

22.  41 

. 2217E*03 

19.  10 

.  2214E*03 

22  49 

. 2214E*03 

19.  17 

. 2210E+05 

22  34 

.  2209E+03 

19  22 

.  221 3E*03 

22  61 

.  2210E+03 

19.  29 

.  2220E*05 

22  68 

.  220BE+03 

19  37 

2219E+03 

22  74 

.  2210E+03 

19  43 

. 2224E*03 

22  81 

.  2209E*05 

19  49 

2227E*03 

22.  86 

221 5E+03 

19.  37 

.  2222E*05 

22.  94 

. 2220E+03 

19  64 

.  2213E+03 

23.  01 

.  2219E*03 

19  70 

. 2217E*03 

23  06 

. 2223E*03 

19  76 

.  2222E+03 

23  13 

. 2223E+03 

19  82 

. 2217E+05 

23  20 

2223E+03 

X  <  cm  ) 


I 


I 


I 


X < cm  ) 

Static  Pr**»oT* 

23  27 

2222E+05 

23.  33 

2224E+05 

23  41 

. 2224E+05 

23  46 

2235E+05 

23  34 

2243E+05 

23  61 

.  2244E+OS 

23.  68 

.  2245E+05 

23  74 

2242E+05 

23  81 

2233E+0S 

23  88 

2243E+03 

23  94 

. 2249E+05 

23  99 

. 2249E+03 

24  08 

. 22S3E+05 

24  13 

.  2250E+05 

24  20 

. 2243E+05 

24  27 

.  2241E+05 

24  33 

.  224BE+05 

24  43 

2261 £+05 

24  30 

2256E+0S 

24  37 

2251E+05 

24  63 

.  2253E+05 

34  71 

2257E+05 

24  78 

. 2237E+05 

24  86 

.  2233E+03 

24  92 

.  2249E+05 

23.  OO 

.  2243E+05 

23  06 

22S0E+05 

23  14 

2254E+OS 

23  22 

2268E+0S 

23.  29 

.  22626+03 

23  36 

22B0E+05 

23  44 

.  2285E+03 

23.  30 

2274E+03 

23  se 

2272E+05 

25  66 

2273E+0S 

23  73 

. 2283E+03 

23  80 

2304E+03 

23  88 

.  2312E+03 

23  93 

231  1E+05 

26  02 

. 2310E+05 

26  09 

2319E+05 

26  16 

.  2326E+03 

26  23 

. 2323E+05 

26  30 

.  2300E+03 

26  42 

. 2302E+03 

26.  43 

. 2299E+03  . 

26.  34 

.  2292E+03 

26  60 

.  22B0E+05 

26.  64 

2279E+03 

26  72 

.  2274E+03 

26  79 

2278E+03 

26  83 

22B4E+03 

26  94 

2290E+03 

27  02 

22B4E+03 

27  08 

22B2E+03 

27  13 

2284 E +05 

27  23 

2290E+03 

27  30 

2296E+05 

27.  36 

.  2304E+0S 

27  44 

. 2314E+03 

27.  31 

.  2318E+05 

Static  Pr«i* 


27  38 

2311E+05 

27  66 

23C5E+05 

27  73 

2299E+0? 

27  79 

2298E+05 

27  87 

2325E+C-3 

27  93 

2301 E+05 

28  02 

2298E+05 

28  08 

2276E+05 

28  13 

2233E-05 

28  23 

2246E+03 

28  30 

. 2236E+03 

28  37 

2263E+0? 

28  48 

2233E+03 

28  33 

. 2238E+03 

26  62 

. 2227E+03 

28  69 

2222E+03 

28  76 

221 1E+05 

28  S3 

2204E +05 

28  91 

.  2199E+03 

28  98 

2200E+03 

29  06 

.  2200E+03 

29  12 

2195E+03 

29  19 

21986+03 

29  27 

2209E+05 

29  33 

2204E+05 

29  41 

21 93E+03 

29  48 

2191E+05 

29  36 

2184E+05 

29  62 

. 2180E+05 

29  70 

217BE+05 

29  77 

.  2176E+05 

29  84 

21766+03 

29  91 

. 2176E+03 

29  98 

.  2171E+0S 

30  06 

21746-05 

30  12 

21786+05 

30  20 

2170E+05 

30  27 

2162E+05 

30  34 

.  2166E+05 

30  40 

2170E+05 

30  48 

. 2177E+05 

30  33 

.  217BE+05 

30  63 

217BE+05 

30  70 

. 2179E+05 

30.  77 

. 2176E+03 

30  84 

. 21 74E+05 

30  91 

. 217ZE+05 

30  98 

.  2173E+05 

31  03 

. 2175E+05 

31  13 

21B2E+05 

31  19 

2194E+05 

31.  23 

. 2202E+05 

31.  32 

2206E+05 

31  40 

.  2218E+05 

31  47 

. 2222E+05 

31  33 

224 1E+05 

31  61 

2262E+05 

31  67 

226BE+05 

31  74 

2273E+05 

31  82 

. 2291E+05 

31  88 

2292E+05 

Static  Prtnurt 


X  (cat.  > 


Static  Prmuri 


Xlcin  ) 

31  94 

32  00 
32.  06 
32.  14 
32  10 
32  29 
32.  30 
32  36 
32  42 
32.  40 
32.  53 
32.  59 
32.  65 
32  71 

32  76 
52  e2 
32.  00 
32.  95 

32.  99 

33  06 

33.  1 1 
33  17 
33  23 
33.  29 
33  35 
33  41 
33.  47 
33  52 
33.  59 
33  64 
33.  70 
33.  75 
33  01 
33  98 

33.  93 

33  99 

34  05 

34  io 

34  16 
34  22 

34  28 
34  34 

34.  40 
34.  46 
•34  52 

34.  58 
34  64 
04.  69 
34  75 
34  81 
34  07 
34  94 

34  99 
'15  05 
25  1  1 

35  17 
35  22 
35  29 
35  34 
35  40 

35.  46 
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2302E+0S 
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.  2273E4-05 
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. 2277E*05 
.  2276E4-03 
.  2276E*05 
. 220OE+O3 
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2283E+05 
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. 2287E+05 
.  2284E«-05 
2285E+05 
.  2284E+05 
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2280E+05 
. 2280E+05 
2282E+05 
2280E-405 
2277E+05 
2280E-405 
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2283E+05 
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2289E+05 
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. 2305E+05 
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37  97 

38  03 
38  09 
38  15 
38  20 

38.  27 
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2376E*03 

44  40 
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2251E-C5 
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2239E+05 

41  29 
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. 2306E+03 

43  73 
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2237E*05 

47  40 
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22386*03 

47  46 

. 23108*05 

31  37 

22416*03 

47  S3 

23206*03 

31  44 

2239E*03 

47  (  60 

23226*03 

31  31 

.  2236E*03 

47  69 
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48  19 
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43  26 
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48  33 

23176*03 
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39 
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48.  46 

23076*03 

32  32 
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23046*03 

32  40 

230 IE *03 

48  60 

23036*03 

32  46 

2296E*03 

48  66 
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2297E*03 

48  72 

23O0E*03 

32  38 
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48  79 
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32  65 
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48  88 
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2276E*03 
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49  Ol 
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49  14 
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33  04 
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49  21 

2237E*0S 

33  11 
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33  17 
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49  33 

. 2236E*03 

33.  23 

.  2281 E*03 

49.  40 

.  2261E*03 

S3  29 

.  2281E*03 

49  45 

2238E*05 

33  37 

. 2280E*03 

49  52 

2231E+03 

S3.  42 

.  2280E*03 

49  58 

. 2246E*03 

33.  49 

. 2279E*03 

49  65 

2242E*03 

33.  33 

. 2279E*03 

49  71 

. 2239E+03 

S3  61 

.  2281E*03 

49  77 

2240E*03 

33  66 

. 2282E+03 

49  84 

. 2241E*03 

33  74 

. 2283E*03 

49  90 

. 2243E*03 

S3  BO 

2284E*03 

49  97 

.  2235E*03 

33.  87 

.  2283E*03 

SO  04 

. 2228E*OS 

33  93 

.  2284E*03 

30  10 

. 2224E*03 

33.  99 

. 2284E*03 

50  IS 

. 2223E*03 

34  03 

.  2285E*03 

50.  22 

. 2227E*03 

34  12 

.  2289E*03 

50  28 

2227E*03 

34  17 

2283E*03 

VJ  OS 

. 2229E*03 

54  24 

.  2284E*0S 

SO  41 

2229E+03 

34.  31 

.  2286E*03 

50  48 

2223E*03 

34  38 

.  2286E*05 

50.  54 

. 2220£*03 

34.  44 

2284E*03 

30  60 

2220E*03 

34  SO 

2284E*03 

“0  o7 

2223E  «-03 

34  36 

2284E*03 

SO  73 

2226E*03 

34  63 

22S6E*03 

50  81 

2230E*03 

34  69 

2284E*03 

X  <  cm  > 

Static  Prttti 

t4  ei 

22B5E+09 

94  ee 

2284E+09 

94  44 

2203E+OS 

99  01 

. 2284E+05 

95  07 

22S6E+05 

99  13 

2284E+09 

99  14 

.  2292E+05 

99  32 

.  2299E+05 

99  32 

2300E  *05 

99  34 

231 OE+OS 

99  49 

.  23138+09 

99  91 

.  2316E+05 

59  58 

2329E+09 

99  64 

2331E+09 

95  70 

2340E+09 

55  77 

. 2392E+09 

95  82 

2360E+0S 

95  84 

.  2368E+05 

59  46 

2373E+05 

56  02 

237BE+05 

96  06 

. 2382E+09 

56  15 

2389E+05 

56  21 

. 2385E+05 

56  27 

.  23BBE+05 

96  33 

. 2340E+05 

96  40 

. 2346E+05 

96  47 

.  2391E+09 

96  93 

.  2387E+05 

96  60 

.  2383E+09 

56  66 

. 2379E+05 

56  73 
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96  77 

. 2372E+09 

96  85 
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96  40 

2363E+09 

56  47 

.  2348E+05 

57.  03 

.  2393E+05 

97  12 

.  2344E+05 

97  14 

.  2341E+05 

97  23 

.  2341E+05 

97  30 

2336E+09 

57  38 

. 2346E+05 

57  42 

.  2334E+05 

97  52 
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97  57 

2328E+05 

97  64 
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57  69 

2330E+05 

97  74 

2326E+09 
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97  84 

2334E+05 
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.  2333E+05 

58  00 

2334E+05 

98  08 

2338E+09 

58  15 

2336E+05 
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2341E+09 

58  28 

2343E+09 

54.  33 

2346E+09 

98  40 

.  2359E+05 

58  46 

.  2395E+05 

98  93 

2346E+09 

98  60 

2346E+05 

58  66 

2353E+05 

X  <  c  #•  > 

Static  Pr*M 

98  75 

235  ’*E*05 

98  80 

23fc2E*05 

5e  90 

2357r*C5 

58  93 

2355E-0I 

54  01 

235 of- 

99  08 

2346£*05 

59  15 
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59  23 

2345E+05 

59  31 
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54  36 

2334fc-<'>5 

59  49 

2324E+05 

99  91 

2316E+05 
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2316E+05 

99  69 
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99  73 

. 2310E*05 

59  81 

231 lE+05 

59  86 
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. 2299E*05 

60  00 

2294E+05 
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2295E+05 

..60  14 

22948*05 
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60  29 
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60  35 
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60  43 

2275E-05 

60  48 

2271E+05 
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60  71 
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2240E-05 
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61  26 
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61  30 
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62  48 

2392E+05 
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2393E+05 

62  64 

2390E+05 

62  69 

2388E+05 

62  76 

2380E+05 

62  84 

2373F*05 

62  91 

236®E*05 

X  (  cm  ) 
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63 

06 

2371E+03 
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1 1 

2373E+05 

63 

18 
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63 

25 
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38 
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63 
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63 
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63 

71 
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63 
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63 

34 
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63 

91 
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63 

97 
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04 
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=4 
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64 

18 
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64 

24 
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64 

30 

2293E*C3 

o* 

36 
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64 

44 
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56 
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-a 
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:4 
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_4 
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64 

98 
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64 

96 

22646*03 

63 

01 
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m. 

« *i 
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66 
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23 
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65 

34 
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42 
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6? 

47 
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65 

54 
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65 

60 
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"■O 

66 

22916*05 

65 

72 
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80 
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65 

86 
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63 

92 
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65 

93 
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66 

06 
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VERTICAL  PITOT  PRESSURE  PROFILE  AT  A  STREAMUISE  STATION 
r»  i6  Si  c» 
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. 2190E+03 

470  3 

1  69  V. 

3731 

U63E-06 

2I90E+03 

473  2 

1  =2; 

6704 

. 124SE*06 

2190E+03 

483  6 

2  004 

7666 

.  X333E*06 

2I90E+03 

497  2 

2  1C  1 

.  8666 

. 1395E+06 

. 2190E+05 

SC  1  2 

2  i:<6 

9616 

1490E+0e 

2190E-03 

31 C  1 

2  2-4 

1  036 

. 1346E+06 

. 2190E+03 

314  fc 

2  2t“ 

J  130 

1600E+06 

2190E+03 

319  o 

2  2=  = 

1  246 

1671E*06 

2190E+03 

324  7 

2  33t 

1  343 

. 1743E+06 

.  2190E+03 

330  2 

2  412 

1  442 

1B33E+06 

2190E«-03 

336  1 

2  4-6 

1  639 

. 1962E*06 

. 2190E+03 

344  0 

£  5eS 

1  823 

2061 E*06 

2190E+03 

34E  9 

2  632 

2  018 

. 2169E+06 

.  2190E+03 

334  2 

2.  703 

2.  212 

.  2268E-06 

2190E+03 

539  1 

2  774 

2  411 

. 2367E+06 

. 21 90E*03 

362  7 

2  830 

2  399 

. 2403E+06 

. 2190E+0S 

364  0 

2  859 

2  791 

. 2422E+06 

. 2190E+0S 

364  5 

2  Bnf. 

2  98? 

.  2446E-I-06 

. 2190E*03 

365  6 

2  8‘  : 

3  177 

2440E+06 

2190E+03 

365  fc 

2  860 

3  370 

. 2463E+06 

. 2190E+0S 

366  6 

2  894 

3  371 

. 2447E+06 

2190E+03 

363  6 

2  ss: 

3  764 

. 2473E*06 

.  2190E*0S 

367  1 

2  °vi 

3  932 

2464E-06 

. 2190E+03 

366  6 

2  e'4 

4  133 

.  2469E • 06 

2190E *03 

366  6 

2  99  4 

4  342 

2463E+06 

. 2190E+05 

366  6 

2  eci 

4  334 

. 24  7E*06 

. 2190E+03 

363  6 

2  860 

4  734 

2473E+06 

. 219OE+03 

367  1 

2  901 

4  930 

2446E*06 

. 2190E+03 

365  6 

2  66  . 

3  119 

. 2460E*06 

2190E+03 

366  6 

2  8  =  4 

rg  i  c  rj  rj  rg  r: 


VERTICAL  PITOT  PRESSURE  PROFILE  AT  A  STREAMWISE  STATION 
24  13  c.s 
l-  O  0  cm 


Average  Stagnation  Prmure"  6904E*O6 
Average  Stagflation  Temperature*  290  7 


V'  cm  ) 

Pitot  Preuure 

Static  Pressure 

U(m/t) 

M 

ea90E-02 

5316E*05 

.  2209Er05 

342  3 

1  200 

6I08E-01 

7052Er05 

?209E«-03 

392  0 

1  441 

1027 

8050Er05 

2209E*05 

414  1 

1  361 

1 4*o 

B704Er05 

2209E+05 

426  4 

1  632 

:  >  33 

9222E-05 

. 2209E*05 

433  7 

1  o89 

-362 

9687E«-05 

2209E*05 

443  6 

1  738 

2876 

1034E>06 

2209Er05 

433  6 

1  304 

3360 

106lEr06 

. 2209E*05 

457  7 

1.  832 

3814 

1 1 00E*06 

2209Er05 

462.  8 

1  867 

•••?'« 

1147Er06 

. 2209E*0S 

468  7 

1  910 

iZO*r 

1 177Er06 

2209E+03 

472  5 

1  938 

3728 

1234E-06 

2209E+0S 

479  0 

1  9B8 

6715 

1294E‘06 

. 22C9E*05 

486  1 

2.  044 

7732 

1386E*06 

. 2209E+03 

493  6 

2.  122 

.  ZL-ib 

143eE*06 

. 2209E*03 

500  4 

2  163 

9586 

1 522E-06 

2209E-05 

507  4 

2  228 

1 .  053 

1583E*06 

.  2209E*05 

512  7 

2  278 

1  1T.3 

165~Er06 

2209E+05 

518  4 

2  333 

l  249 

1 /55E*06 

.  2209E+05 

525  4 

2  403 

’  7  iO 

1816E-06 

. 2209E*05 

529.  9 

2  433 

1  442 

ia.-9E*06 

2209E*05 

533  6 

2.  497 

i  631 

. 203SE*06 

. 2209E*05 

545  7 

2  639 

1  324 

.  22G6E  *-06 

2209E+05 

551  4 

2  717 

2  020 

. 2331 Er 06 

. 2209E*05 

556  8 

2  793 

21a 

245CE*06 

. 2209E+09 

561  9 

2  873 

2  407 

2473E+06 

.  2209E+05 

562  1 

2  887 

2.  596 

. 2497E*06 

2209E+03 

562  3 

2  901 

2  787 

2514E-06 

.  2209E+05 

562  8 

2  908 

2.  983 

.  2343E+06 

2209E*05 

564  3 

2  929 

3.  130 

2346E+06 

. 2209E+05 

564  3 

2  929 

3.  374 

2344E+06 

. 2209Er05 

564  3 

2  929 

3  567 

2344E-06 

.  2209E*05 

564  3 

2  929 

3.  ■’f-a 

.  2341E*C6 

2209E*05 

564  3 

2.  929 

3  434 

2539Er06 

. 2209E*05 

563  8 

2  922 

4  1 52 

.  2363E+06 

. 2209E*03 

565  3 

2  943 

4  343 

.  25B5E*06 

. 2209E*05 

565  B 

2  951 

».  5.73 

.  2393E+06 

. 2209E+03 

566  3 

2.  938 

4  /27 

2374E*06 

. 2209E*03 

565  3 

2  943 

4  724 

2556Er06 

. 2209E*05 

564  8 

2.  936 

5  122 

253 1 E» 06 

2209E*03 

363  8 

2  922 

VERTICAL  PITOT  PRESSURE  PROFILE  AT  A  STREAMWISE  STATION 
X*  31  75  cm 
2*  0  0  cm 


Avtngt  Stagnation  Prruuri*  6913E*06 
Average  Stagnation  Tampanaturam  250  0 


Y<  cm  ) 

Pitot  Praiiura 

U  t  m  /  i  ) 

m 

839u£-02 

.  5145E+05 

2l9lE*05 

337  1 

1  1-9 

64B6E-01 

7326E+05 

2191E+0S 

399  5 

1  a-  :■ 

1076 

B09BE*05 

21 9 1E*05 

415  4 

1  5-2 

.  1579 

8795E+-05 

. 2191E+05 

426  7 

1  c ; 

.  2074 

9302E*05 

2191E+05 

437  9 

3  7r>. 

2494 

9705E*05 

2l 93  E*05 

444  0 

1  7_* 

3033 

1021E*06 

21  91E*05 

4*  1  g 

1  ' 

3476 

.  1051E*06 

. 2191E+05 

456  0 

1  S 7 

3965 

.  1106E*0o 

21 91E*05 

4C4  1 

1  8=: 

4487 

.  1 1 52E+06 

. 21918*05 

470  0 

1  924 

4926 

.  1 1 79E*06 

. 2191E+05 

472  e 

1  94'. 

5894 

.  1243E  +  06 

2191E+05 

4g  1  1 

2  0’>° 

6897 

.  130?E*G6 

. 2193E+05 

483  1 

2  jft 

7046 

13828*06 

. 2191E+05 

495  7 

2  12- 

e780 

1429E*06 

. 2191E*05 

4=9  6 

2  i  -  * 

9737 

.  14798*06 

. 2191E*05 

504  r 

2  2'*- 

1  073 

1544E+06 

. 2191E+05 

5  5 

2  2*7 

1  169 

1653E*06 

. 21 91 E*05 

5 1  S  4 

2  34.7 

1  263 

.  1685E+06 

. 2I91E+05 

52C  2 

2  3c  2 

1  35c 

.  1 794E+06 

. 2191E*05 

522  4 

2  44  = 

1  452 

. 1836E*06 

.  2191E*05 

530  8 

2  4-; 

1  647 

.  2004E*06 

2191E+05 

54  3  1 

•■t  «  C  ”7 
m.  w  ■*  / 

1  844 

. 21 60E*06 

. 2191E+05 

549  4 

2  702 

2  036 

2252E+06 

. 2191E+05 

552  2 

7  75  = 

2  226 

.  2379E+06 

. 2191E+05 

55c  4 

2  E"’’ 

2  43(7 

.  2446E*06 

. 21 91 E*05 

560  E 

2  so: 

2  609 

2504E+06 

. 2191E*05 

562  5 

2  9:5 

2  808 

. 2523E+06 

. 2191E+05 

56  3  6 

2  92  = 

3  001 

. 252BE+06 

. 2I91E+05 

562  6 

2  92= 

3  189 

. 2546E+06 

21 91E*05 

564  6 

2  94  2- 

3  382 

. 2529E*06 

. 2191E+05 

56?  6 

2  92  * 

3  583 

.  2534E*06 

. 21 9 i£*05 

564  ) 

2  9  36 

3  779 

2539E+06 

.  2J91E+05 

564  1 

2  ®3t 

3  971 

.  25248*06 

2l91E*05 

5o3  6 

2  =2= 

4  162 

. 2532E*06 

.  2191E+05 

564  1 

2  9Jc 

4  354 

. 2521E*06 

21 9iE*0S 

562  6 

2  92= 

4  551 

2519E+06 

2191E+05 

563  6 

2  92= 

4  748 

. 252SE*06 

. 2191E+05 

562  6 

2  92= 

4  940 

25228*06 

.  19IE+0S 

56?  6 

2  92= 

5  129 

2 542E*06 

2191E+05 

5  64  J 

2  93c 

'v  Eft  T I  CAL  PITOT  PRESSURE  PROFILE  AT  A  STREAHWISE  STATION 
<*  39  3?  cm 
;»0  0  cm 


Avaraga  Stagnation  Prttiurt*  . 6903E+06 
Avaraga  Stagnation  Taaiparatura*  291  1 


V  ( cm  ) 

Static  Pra*»ur* 

U(m/«  > 

239CE-02 

.  9269E*09 

2267E*09 

336  3 

5435E-01 

7151E*09 

2267E*09 

390  3 

9701E-01 

.  7869E+09 

2267E*09 

406  2 

1473 

8401E-*-O9 

2267E+09 

418  7 

1916 

B973E*09 

2267E+09 

427  8 

2434 

971GE-09 

. 2267E+09 

439  9 

*  'j'&Q 

101 9E*06 

2267E*09 

447  1 

2279 

1049E+06 

2267E+09 

491  6 

2909 

1077E~O6 

2267E+09 

499.  8 

4332 

1 084E*06 

. 2267E+09 

496  7 

4  370 

1127E+06 

. 2267E-09 

462  a 

;ai9 

1224E*06 

2267E*09 

474  9 

i’23 

131 9E*06 

2267E+09 

489  6 

' 7 1  3 

1347E+06 

2267£-»09 

488  O 

3651 

1412E*06 

2267E+09 

494  7 

9657 

1 900E*06 

2267E*09 

902  1 

:  066 

1947E-M36 

. 2267E+09 

906  6 

:  : *4 

1 97 1 E*06 

.  2267E+09 

907  9 

1  240 

1633E*06 

.  2267E+09 

916  9 

1  J44 

1740E-KJ6 

2267E+09 

921.  0 

:  447 

1789E+06 

.  2267E+09 

924  3 

:  635 

1870E*06 

2267E*09 

929  2 

:  327 

2087E*06 

.  2267E+09 

942  2 

v  0 22 

2243E-M26 

2267E+09 

949  9 

*■  2^2 

. 2342E*06 

2267E*09 

994  2 

e  4 :e 

. 24O9E+06 

.  2267E*09 

996  7 

i  597 

2467E*06 

.  2267E+05 

998  6 

791 

2490E-06 

.  2267E+09 

999  6 

2  997 

2906E+06 

.  2267E+09 

960  2 

7  1 96 

2499E*06 

.  2267E+09 

999  6 

J  373 

2912E+06 

2267E+09 

960  7 

3  5-4 

2903E*06 

.  2267E709 

960  2 

7  767 

. 2920E*06 

.  2267E+09 

960.  7 

7  956 

2944E*06 

2267E+09 

961  7 

4  155 

2939E»06 

2267E*09 

961  7 

4  242 

2996E+06 

2267E*05 

962  2 

4  540 

2999E+06 

2267E+09 

962  S 

4  731 

2942E-06 

2267E*09 

961  7 

'•  929 

2944E+06 

2267E-09 

961  7 

123 

2938E  *06 

2267E*09 

961  7 

H 

1  172 
1  430 
1  919 
1  986 
1  639 
1  713 
1  799 
1  789 
1  818 

1.  829 
1  867 

1  992 

2  037 

2.  098 
2.  119 
2.  179 
2.  221 
2.  239 
2.  320 
2  363 
2  398 
2  499 
2.  604 
2  703 
2  766 
2  809 
2.  844 
2  899 
2.  866 
2.  899 
2  873 
2  866 
2  873 
2  887 
2  887 
2  894 
2  901 
2  887 
2  887 
2  387 
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MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  OF  STANDARDS- 1963-A 


VERTICAL  PITOT  PRESSURE  PROFILE  AT  A  STREAMUISE  STATION 
X*  40.  18  CM 

Z»  0.  0  cm 


Avoragt  Stagnation  Prttiurt*  .  6B97E+06 
Avvrag*  Stagnation  Tamparaturaa  232  3 


Y (cm  ) 

Pitot  Prtttun 

Static  Pra**ur* 

U  ( ■•/  i  ) 

n 

.  8890E-02 

.  9291E+09 

.  2276E4-09 

335  5 

1  165 

. SB06E-01 

-  71 1 1E+05 

2276E4-09 

3B7  e 

1  422 

.  1019 

.  8078E4-09 

.  2276E4-09 

411  0 

1  536 

.  1979 

.  879BE4-09 

.  2276E4-09 

42 A  6 

1  6!< 

.  1999 

. 9410E*09 

.  22 7 bE  4-0  9 

435  3 

1 .  6,S 

2494 

.  9923E+05 

r  .  2276E4-09 

443  C 

1  721 

.  2971 

.  9964E*05 

. 2276E+09 

444.  £ 

1  7^v 

.  3409 

.  1094E4G6 

.  2276E4-09 

452  7 

1  tF- 

3977 

.  1089E406 

.  2276E4-09 

456  e 

1  B’S 

.  4423 

.  112BE+06 

.  2276E+05 

463  0 

1 .  BsC 

.  4937 

.  1 1 73E*06 

.  2276E4-09 

469  0 

i  9c: 

.  9879 

.  1292E«06 

.  2276E4-09 

47E  f 

1  974 

6817 

.  1307E»06 

2276E4-09 

484  e 

2  023 

.  7791 

.  1349E4-06 

. 2276E+09 

48E  2 

2  OM 

.  8716 

.  1444E«-06 

2276E4-09 

49E  5 

2  I2e 

9707 

.  1481E*06 

2276E4-09 

901  6 

2.  165 

1  070 

.  1949E4-06 

. 2276E+09 

507  C 

2  ?>4 

1.  166 

.  1367E-06 

.  2276E4-09 

508  3 

2  22= 

1. 297 

.  16B0E+06 

.  2276E4-09 

917  4 

2.  312 

1  349 

.  171 4E*06 

.  2276E4-09 

920  0 

2  347 

1.  491 

.  1809E«06 

.  2276E4-09 

926  2 

2  405 

1  643 

.  1979E-*06 

.  2276E4-09 

937  2 

2  52= 

1.  837 

.  2088E*06 

.  22 7 6E  4-0  9 

542  e 

2.  5'-7 

2  027 

.  2146E*06 

.  2276E*09 

545  8 

2  62c 

2  222 

.  2346E+06 

227 6E 405 

955  9 

2  76 * 

2.  420 

.  2410E+06 

.  2276E405 

557  5 

2  802 

2  609 

.  2442E*06 

.  2276E409 

958  3 

2  822 

2.  602 

2923E406 

.  227 6E 409 

962  0 

2  S"2 

2.  993 

. 2926E*06 

.  227 6E 409 

962  0 

2  872 

3.  189 

. 2919E+06 

. 2276E405 

561  9 

2.  e.-i 

3  389 

. 2946E*06 

.  2276E409 

563  1 

2  err 

3  979 

. 2934E*06 

. 2276E409 

562  6 

2  BBC 

3.  772 

.  2939E4-06 

.  2276C409 

562  6 

2  BSC 

3.  960 

.  2931E  +  06 

.  2276E409 

562  0 

2  e"’2 

4  199 

.  29O9E+06 

. 2276E409 

561  0 

2  8*  ‘ 

4  393 

. 2929E+06 

.  2276E409 

562  0 

2  e72 

4.  948 

.  294  3E 4-06 

.  2276E409 

962  6 

2  at* 

4  738 

.  2941E4-06 

.  227 6E 409 

962  6 

2  esc 

4  931 

. 2953E+06 

.  2276E409 

563  1 

2  8P~ 

9  126 

.  2999E406 

■  2276E409 

563  6 

2  8«4 

VERTICAL  PITOT  PRESSURE  PROFILE  AT  A  STREAMUISC  STATION 
(•  46.  99  cm. 

Z»  0.  O  cm. 


Baraga  Stagnation  Prtnurt*  6801 E* 06 


-varaga  Stagnation  Tanpa*atura» 

292  8 

V(C*.  > 

Static  Preaaura 

U(m/a  > 

H 

33801-02 

4621E*09 

.  92939*03 

313.  6 

1.  0*8 

52281-01 

. 6690E»09 

92339*03 

381  I 

1.  377 

90=»6E-01 

. 7922E-09 

33339 *09 

401  0 

1  480 

1242 

.  S399E-K39 

2292E*09 

418  3 

1  979 

:*24 

. 8/881*09 

2292E*09 

42*  4 

1  *29 

2010 

9183E-09 

2292E*09 

434  0 

1  4*7 

244? 

.  89021*09 

. 2292C*09 

438  t 

1  *88 

2916 

.  8890E*09 

2292E*09 

444  7 

l  739 

2156 

1090E-06 

33339*03 

494  3 

1  787 

*5*9 

10831*06 

2292E*09 

460  4 

t  838 

2?  2\ 

1082E*06 

. 2292E*09 

4*0  9 

l  838 

U46C*06 

2292E*09 

467  9 

1  888 

5524 

12091*06 

2292E*09 

479  2 

1  849 

6250 

12461*06 

2292E*09 

478  8 

1  881 

*039 

12861*06 

2292E*09 

489  2 

2  023 

**59 

1 349E*06 

2292E*09 

480  9 

2  06* 

*597 

1403E*C6 

2292C*09 

48*  4 

a  ns 

^340 

14661*06 

2292E*09 

902  2 

a  1*9 

1  006 

19361*06 

33339*03 

908  * 

2  221 

.  OPS 

1949E*06 

33339*03 

908  1 

2  228 

•  169 

16411*06 

.  33339*03 

916  7 

2.  288 

:  316 

1704E*06 

23339*03 

921  4 

2  348 

■  4 ’3 

18?ll*C6 

2292E*09 

931  6 

2  499 

-27 

18181*06 

2292E*09 

939  1 

3  487 

1  777 

18881*04 

22921*09 

938  8 

2  994 

*»29 

20731*06 

33339*03 

943  * 

3  *04 

'  099 

22201*06 

2292E*09 

991  9 

3  703 

Z  247 

32691*06 

22921*09 

993  1 

3  731 

:  4oo 

23 4 BE -0* 

2292E*09 

99*  4 

2  781 

2  948 

23701 *04 

22921*09 

99*  8 

3  789 

?  yss 

24381-04 

22921*09 

998  8 

3  837 

2  060 

34381*06 

22921*09 

998  8 

2  837 

3  017 
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FIGURE  1 


The  Mechanics  of  streak  breakup;  adapted  from  Kline 
et.  al  (1967) 


FIGURE  2. 


Side  view  of  the  Interactions  between  bursting  flow 
eodu1es(eode1  Interpretation)  (Offen  and  Kline,  1975) 


Schematic  showing  the  development  and  Interaction  of 
hairpin  vortices  created  by  the  three-dimensional 
separation  of  the  flow  over  a  hemisphere  (Acarlar 
and  Smith,  1984) 
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c)  vortoa  dovolopoant:  ossification  and  cone ant  ration 
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d)  vortoa  aj action,  stretching  and  Intoraction 


FIGURE  4. 


Illustration  of  the  breakdown  and  formation  of 

hairpin  vortices  during  the  bursting  process.  Low* 

spaad  streak  regions  are  Indicated  by  shading  (Smith ,1984) 


FIGURE  5. 


Effect  of  Reynolds  nueber  on  eddies:  (a)  very  low  Re 
(vortex  loops),  (b)  lownaoderate  Re ( horseshoes ) , 

(c)  laode  rate-high  Re(halrplns)  (Head  and 
Bandyopadhyay ,  1981) 


FIGURE  6.  20°  Interface  caused  by  45°  structures 

(Head  and  Bandyopadhyay ,  1981) 


FIGURE  7. 


Random  array  of  eddies  at  a  constant  angle  The 
sampling  volume  for  the  p.d.f.  Is  ABCO  (Perry 
and  Chong,  1982) 


Stk  Wwntfcy 


FIGURE  8. 


Symbolic  representation  of  a  discrete  system  of 
geometrically  similar  hierarchies  (Perry  and 
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Sketch  of  the  Kullte  pressure  transducer  plug 
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FIGURE  10  b  Arrangement  of  hot  wire  probes  and  wall  pressure  transducers 
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FIGURE  13. 


Flowchart  of  the  hot-wire  data  acquisition  system 
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Normalized  Probability  Density 
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FIGURE  15.  The  mass-flow  probability  density  function  at 

several  points  In  the  boundary  layer 
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FIGURE  19. 
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Typical  naan  velocity  profile 
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FIGURE  21. 
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Spenwlse  Pitot  pressure  distributions  «t  *  • 
18.5  In.,  y  »  0.15  end  0.4  In. 
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Typical  r.m.s.  mass-flow  fluctuation  level  through 
the  boundary  layer 
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FIGURE  24. 
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FIGURE  26. 
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Broad  band  convection  velocities,  deduced  from 
wall-pressure  signals 
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FIGURE  28. 
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The  computed  wall -pressure  coherence  functi 
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FIGURE  31 


The  computed  wall -pressure  time-delay  function. 
See  Figure  29  for  legend 
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FIGURE  32. 


The  computed  wall  pressure  convection  speed  function. 
See  Figure  29  for  legend 


FIGURE  34a. 


Sketch  of  a  typical  large-scale  structure.  Inclined 
at  an  angle  0,  passing  a  double  wire  probe 


FIGURE  34b.  Diagram  used  to  determine  "structure  angle  " 
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FIGURE  35. 


Computed  large-scale  structure  angle  throughout  the 
boundary  layer  for  two  different  wire  separation 
distances 
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Coherence  function  (through-dut  the  boundary  laydr) 
for  two  fluctuating  mass-flow  signals  separated  by 
0.096. 


FIGURE  37a. 


Space-time  correlation  between  a  wall-pressure 
signal  and  a  mass-flow  signal.  Streamwlse 
separation  of  0.456 
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FIGURE  37c.  The  same  as  Figure  37a;  streamwise  separation  of 

0.096,  spanwise  separation  of  0.116 
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FIGURE  38.  Instantaneous  wall-pressure  signal  and  Its  corresponding 

short-time  variance 
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FIGURE  41. 
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Effect  of  the  threshold  level  on  the  average  posltl 
pressure  event 
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Effect  of  the  short-time  variance  period  on  the  average 
positive  pressure  event 
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FIGURE  43.  "Optimized"  average  positive  pressure  event 
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FIGURE  44.  "Optimized"  avenge  negative  pressure  event 
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FIGURE  46.  Average  negative  pressure  zerocrossing 


FIGURE  47.  Instantaneous  mass-flow  signal  and  Its  corresponding 

short-time  variance 
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FIGURE  49. 


Frequency  of  mass-flow  events  in  the  boundary  layer 
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FIGURE  51. 
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Effect  of  the  threshold  level  on  the  frequency  of  mass 
flow  events 


FIGURE  53. 


"Optimized"  average  negative  mass-flow  events  at  three 
points  through  the  boundary  layer 


FIGURE  54.  Three  simultaneously  measured.  Instantaneous 

mass-flow  signals  with  their  corresponding  short 
time  variances 
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FIGURE  57. 
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Average  positive  pressure  events  conditioned 
upon  mass-flow  events  from  different  locations 
In  the  boundary  layer.  The  detection  wire  was 
located  at:  A)  y/6  ■  0.03;  B)  y/6  •  0.13; 

C)  y/6  »  0.22;  D)  y/6  »  0.30 


Pressure  transducer  located 
the  hot-wire  probe 
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FIGURE  57  c 


Pressure  transducer  located  0.456  upstream  of 
the  hot-wire  probe 


TIME  MS6C 


FIGURE  57  d  Average  positive  mass-flow  events  from  which  pressure 

events  were  detected 
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FIGURE  58. 


Average  positive  pressure  events  conditioned  upon 
mass-flow  events  from  different  locations  In  the 
boundary  layer.  The  detection  wire  was  located  at- 
J  “  0.03;  F)  y/6  -  0.12;  G)  y/6-  0.21; 

Pressure  transducer  located  1.416  upstream  of  the 
hot-wire  probe 
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Average  positive  mass-flow  events  from  which 
pressure  events  were  detected 


Average  positive  pressure  events  conditioned  upon 
mass-flow  events  from  different  locations  in  the 
boundary  layer.  The  detection  wire  was  located  at 
I)  y/6  -  0.03;  J)  y/6  -  0.08;  K)  y/6  -  0.15; 
L)  y/6  -  0.20 

(a)  Pressure  transducer  located  0.096  upstream  and 
0.116  laterally  from  the  hot-wire  probe 


Average  positive  pressure  events  conditioned  upon 
mass-flow  events  from  different  locations  in  the 
boundary  layer.  The  detection  wire  was  located 
at:  M)  y/6  =  0.03;  N)  y/6  •  0.08;  0)  y/6  «  0.16; 
P)  y/6  *  0.23 

(a)  Pressure  transducer  located  0.096  upstream  and 
0.235  laterally  from  the  hot-wire  probe 
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Average  positive  mass-flow  events  from  which 
pressure  events  were  detected 


FIGURE  62.  Ensemble-average  time  histories  of  the  conditionally 

sampled  streamwise  velocity.  Detection  is  based  on 
smoothed  rectified  high  frequency  u-component  at 
each  y-value  (Thomas  and  Bull,  1983) 
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FIGURE  63. 


Positive  VITA  conditional  averages  of  the  streamwise 
velocity  at  y+*  13  (for  thyee  different  threshold 
levels),  normalized  by  (Ku  rms)^  (Alfredsson  and 
Johansson,  1982) 


FIGURE  64. 


Positive  VITA  conditional  average  of  streamwise 
velocity  at  y+  «  15  with  various  threshold  levels 
(Blackwelder  and  Kaplan,  1976) 
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TABLE  2. 


M. 

2.87 

Ra./m 

6.5  »107 

6 

28  mm 

6* 

6.2  mm 

U. 

565  "Vs 

<pu>. 

479  k9/m* 

Rvall 

3.34  p.s.i. 

Cf 

.00114 

Flow  conditions  on  the  tunnel  centerline  at 
x  ■  18.5  Inches 


Threshold 

r 

No.  of  Events 

0.7 

.18 

90 

0.8 

.14 

69 

1.0 

e 

O 

<0 

44 

Effect  of  the  threshold  level  on  the  number  of 
pressure  events 


Positive  Events  Negative  Events 


\  f  «* 

.03  .19  .13 

.43  .43  .23 

1.00  .29  .10 


TABLE  3. 


Non-dlmenslonallzed  occurence  frequencies  of 
"optimized"  mass-flow  events  in  the  boundary 
layer 
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